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ABSTRACT 


The  wall  temperature  distribution  of  a  combustion  chanter  Is  a 
function  of  the  various  heat  transfer  processes  existing  In  the  cham¬ 
ber  and  annulus.  The  basic  turbulent  conservation  equations  of  mass, 
momentum,  species,  and  energy  are  developed  In  an  effort  to  provide  an 
analytical  rather  than  empirical  method  of  determining  the  temperature 
distribution.  However,  the  inclusion  of  radiation  energy  exchange 
using  the  radiative  transport  theory  and  the  fact  that  the  chamber 
flow  Is  not  one  dimensional  makes  a  closed  form  solution  to  the  pro¬ 
blem  mathematically  Impossible. 

The  literature  written  prior  to  1966  Is  reviewed  In  an  effort  to 
provide  approximate  analytical  methods  for  determining  the  effects  of 
the  various  modes  of  heat  transfer  on  the  wall  temperature.  Specifi¬ 
cally,  the  review  covers  methods  of  evaluating  radiation  exchange  be¬ 
tween  the  chamber  gases  and  wells,  between  different  areas  of  the 
chamber  wall,  and  between  the  chamber  wall  and  the  annulus  casing. 
Also,  the  methods  available  for  determining  the  effects  of  convective 
heating  and  cooling  In  the  form  of  convective  heat  transfer  In  the  ab¬ 
sence  of  secondary  flow,  film  cooling,  transpiration  cooling,  and  the 

X 

effect  of  penetration  jets  are  reviewed.  Conclusions  are  drawn  as  to 
the  possible  areas  In  which  considerable  advancement  In  calculation 
procedures  have  been  made  since  1966. 
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NOMENCLATURE 

Units 

A 

Area 

ft2 

\ 

Turbulent  viscosity 

lbm/ft  sec 

\ 

Turbulent  thermal  conductivity 

lbm-sec/ft' 

B 

Constant  of  proportionality  from 

Eqn.  (4-61 c) 

C 

Carbon  particle  concentation 

part/ft2 

Cmax 

Maximum  carbon  particle 
concentration  in  chamber 

part/ft3 

C/H 

Carbon  to  hydrogen  ratio  in  the 
fuel  by  mass 

r« 

n 

Constant  In  Eqn.  (4-61  a) 

c.v. 

Control  volume 

C1 

Convective  heat  transfer  in  the 
chamber 

Btu/hr  ft2 

C2 

Convective  heat  transfer  in 
annulus 

Btu/hr  ft2 

cf 

Friction  factor 

CP 

Specific  heat  at  constant  pressure 

Btu/lbm  °K 

c 

Speed  of  light 

ft/ sec 

co 

Speed  of  light  in  a  vacuum 

ft/sec 

D 

Diffusion  coefficient 

ft2/sec 

°p 

Einstein's  (1963)  plug  height 

ft 

“a 

Diameter  of  annulus 

ft 

d* 

Hydraulic  diameter 

ft 

eslot 

Energy  flux  through  a  film 
cooling  slot 

Btu/ft2  sec 

IX 


Uni  ts 


Eb 

Black  body  emissive  power 

Btu/hr  ft^ 

% 

Radiation  flux  vector 

Btu/hr  ft 

F 

Radiation  view  factor 

f' 

Glauert's  (1956)  nondimensional 
velocity 

f(A) 

Function  of  wavelength  and  particle  size 
and  distribution 

Gv 

Radiation  quantity  defined  in  Eqn.  (2-22) 

Btu/hr  ft2 

Gravity  or  body  force  in  ith  direction 

ft/sec2 

Hv 

Radiation  quanity  defined  in  Eqn.  (2-14) 

Btu/hr  ft 

hc 

Convective  heat  transfer  coefficient 

Btu/ft2  sec 

h‘, 

hp 

Defined  as  =  0.0265  (Re)°^S(Pr )0*3 

Btu/ft2  sec 

Planck’s  constant 

Ibm  f t2/sec 

h 

Specific  enthalpy 

Btu/ft3 

h; 

Slot  height 

ft 

>bv 

Black  body  emitted  intensity  at 
frequency  v 

Btu/ft2  hr 
steradians 

!v 

Intensity  of  radiation  at  frequency  v 

Btu/ft2  hr 
steradians 

dIv 

Differential  change  in  I 

Btu/ft2  hr 
steradians 

K 

Conductive  heat  transfer 

Btu/hr  ft2 

K1 

Constant  in  the  flame  emissivity 

Eqn.,  (Eqn.  (3-28) 

K1 

Constant  in  Eqn.  (4-6 lb) 

K1 

Constant  in  Eqn.  (4-64) 

K1 

Constant  in  Eqn.  (4-67) 

kb 

Boltzman's  Constant 

Bt’j/mole  °K 

Ki 

Thermal  conductivity 

Btu/ft  hr  °l 

L 

Width  of  film  cooling  slot 

ft 

Units 

LP 

Einstein's  (1963)  plate  v/idth 

ft 

i 

Mean  beam  length 

ft 

l 

Mean  free  path 

ft 

m 

Nondimensional  mass  flow  rate  through  wall 

“s 

flondimensional  mass  entrainment  from  the 
frees  treani 

m 

Mass  flow  rate 

Ibm/sec 

mOIT 

Rate  of  mass  entrainment  from  mainstream 
into  the  boundary  layer 

lbm/ft2  sec 

ncr 

A  measure  of  the  importance  of  conductive 
to  radiative  heat  transfer  in  a  system 

n 

A  measure  of  the  rate  of  heat  transfer 
and  mass  transfer  to  momentum  transfer 
in  a  system 

Pr 

Prandtl  number 

P 

Pressure 

lbf/in2 

P 

pressure 

lbf/in2 

P 

Time  averaged  pressure 

lbf/in2 

Pv(s-s) 

Probability  function  used  in 

Eqn.  (2-10) 

Q 

Mass  flow  rate 

1 bm/sec 

qC 

Conductive  heat  transfer  flux 

Btu/ft2  sec 

qR 

Radiative  heat  transfer  flux 

Btu/ft2  sec 

qRC 

Combined  radiation  and  conduction  heat 
transfer  flux 

Btu/ft2  sec 

qco 

Convective  heat  transfer  flux 

Btu/ft2  sec 

qri 

Radiation  heat  transfer  per  unit  area 
in  the  itfl  direction 

Btu/ft2  sec 

«; 

Heat  generated  per  unit  volume 

Btu/ft"*  sec 

Ri 

Radiation  flux  from  the  flame  to  the  wall 

Btu/ft2  sec 

R2 

Radiation  flux  from  the  combustor  wall 
to  the  casi ng  wall 

Btu/ft2  sec 

xl 


Re 

r 

r 

max 


\ 

| 


t 


ui 

v 

V 


I- 

j 


i 


x 


Radiation  flux  between  chamber  walls 

Uni  ts 

Btu/ft2 

Reynolds  number 

Local  fuel  to  air  ratio 

Maximum  radius  considered  for 

radiating  hemisphere  above  a  surface 
element 

ft 

Radius  of  emitting  hemisphere  of  gas 

ft 

Radial  length  of  flame  element 

ft 

Radius  of  chamber 

ft 

Stanton  Humber 

Distance 

ft 

Direction 

Direction 

Inlet  temperature  to  the  coK>ustor 

°K 

Rise  In  Temperature  due  to  combustion 

°K 

Temperature 

°K 

Time  averaged  turbulent  temperature 

°K 

Time 

Reference  velocity  for  wall  jet  flow 

sec 

ft/sec 

Velocity  in  1th  direction 

ft/ sec 

Time  averaged  turbulent  velocity 

ft/sec 

Velocity  In  x-di  recti  on 

ft/sec 

Velocity  of  fluctuation 

ft/sec 

Volune 

Velocity  In  y-direct^n 

ftJ 

ft/sec 

Time  averaged  production  rate  of 

species  m  per  uni'  olune 

Tbm/ft^  sec 

Characteristic  length  of  system 

ft 

xll 


Units 

X1 

Direction 

Xfw 

Distance  from  flame  to  wall  element 
for  two  dimensional  radiation  exchange 

ft 

xa 

Distance  from  an  arbitrary  datun  to 
wall  element  on  wall  a 

ft 

xb 

Distance  from  an  arbitrary  da  tun  to 
wall  element  on  wall  b 

ft 

“f 

Elemental  length  of  flame 

ft 

xco? 

Mole  fraction  of  CO2 

o 

CM 

Mole  fraction  of  I^O 

Mass  fraction  of  species  m 

0-zE) 

Measure  of  magnitude  of  free-mixing- layer 
conponent  of  the  velocity  profile 

a 

Thermal  dlffusivity  k/p  c 

ft2/sec 

a' 

Angle  between  s  and  s' 

Degrees 

“w 

Absorptivity  of  the  wall 

af 

Absorptivity  of  the  flame 

aWJ 

A  measure  of  the  extent  of  frictional 
effects  on  the  flow  for  the  velocity 
profile  beyond  the  maximum  velocity 

6 

Coefficient  of  thermal  expansion 

°K-1 

Monochromatic  extinction  coefficient 

ft"1 

r 

i 

Transmittance  of  thermal  radiation 

1 

r 

mean 

Transmittance  of  thermal  radiation 
evaluated  at  mean  chamber  radius 

6 

Boundary  layer  thickness 

ft 

abl 

Boundary  layer  thickness 

ft 

6xa 

Wall  element  length  on  wall  a 

ft 

6x 

Ah 

Wall  element  length  on  wall  b 

f 

ft 

D 

£ 

Emissivity  of  radiating! matter 

Units 


Nonluminous  component  of  emissivity 
Turbulent  eddy  dlffuslvlty 
Nondimensional  distance  from  the  wall 
Monochromatic  emission  of  matter 
Monochromatic  coefficient  of  absorption 
Luminosity  factor 

Function  which  depends  on  the  distance 
from  maximum  velocity  to  the  wall 
divided  by  distance  from  maximum  to 
free  stream  velocity  for  the  wall  jet 

Wavelength 

Viscosi ty 

Frequency 

Densl  ty 

Stefan-Boltzman  Constant 

Monochromatic  coefficient  of 
scatteri ng 

Shear  stress 

Optical  Thickness 

Time  averaged  Rayleigh  dissipation  term 

Angle  of  inclination  from  normal  of  a 
surface  element 

Conserved  property 

A  function  of  r  and  T 

A  function  of  p 

A  function  of  T*.  or  T,  in  emissivity  or 
absorptivity  respectively 

! 

s Solid  angle  along  s 
Solid  angle  along  s' 


lbm  ft/sec 
Btu/hr  ft2 


y 

lbm  ft/sec 
sec"1 
lbm/ft3 
Btu/sec  ft2 


lbf/f r 


sec 


Radians 


Steradians 

Steradians 


Subscripts 


A 

Annulus 

dA] 

Elemental  area 

dA2 

Elemental  area 

CA 

Casing 

f 

Properties  evalucated  at  the  film 
temperature,  (T^  +  T^)/2 

f 

Flame,  for  radiation  quantities 

f-w 

Flame  to  wall 
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Gas 

max 

Maximum  value 

ms 

Mainstream 
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s 
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At  frequency  v 
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CHAPTER  I 


The  proposed  use  of  the  gas  turbine  power  plant  for  vehicular  pro¬ 
pulsion  has  focused  interest  on  the  temperature  distribution  along  the 
combustor  chamber  walls.  Unusually  high  wall  temperatures  are  expected 
In  certain  applications  due  to  the  necessity  of  high  combustor  inlet 
temperatures  to  increase  turbine  efficiency  and  decrease  specific  fuel 
consumption;  this  Increase  in  inlet  temperature  is  the  result  of  use  of 
a  regenerative  cycle  in  vehicular  gas  turbine  engines. 

Due  to  these  unusually  high  temperatures,  an  analytical  approach 
to  the  determination  of  the  wall  temperature  distribution  is  desirable 
for  two  reasons.  Firstly,  the  presently  available  empirical  data  used 
for  predicting  wall  temperature  distributions  will  not  be  applicable  in 
the  high  temperature  ranges.  Thus,  it  would  be  a  significant  accomplish¬ 
ment  If  theoretical  results  could  be  extended  into  the  high  temperature 
region.  Secondly,  as  the  wall  temperatures  increase,  the  necessity  for 
effectively  using  the  available  coolant  gases  becomes  more  important. 
Therefore,  the  current  study  Is  concerned  with  possible  analytical 
analyses  of  heat  transfer  in  gas  turbine  combustion  chambers. 

This  report  is  the  first  portion  of  a  literature  review  which  will 
provide  an  understanding  of  the  empirical  and  theoretical  tools  available 
for  the  analysis.  One  such  analysis  has  already  been  completed  by 
Northern  Research  and  Engineering  Corporation  (Anon.,  1968)  with  limited 
useful  results.  Since  this  analysis  was  completed  in  1966  it  forms  a 
suitable  division  for  the  present  literature  review.  The  first  segment 
of  this  review  concerns  those  articles  directly  used  in  the  NR EC  analysis 
and  other  significant  articles  written  before  1966,  The  second  part  of 
the  review  will  deal  with  those  works  published  since  I960.  The  latter 
will  hopefully  provide  the  background  material  necessary  for  understanding 
and  judging  those  advancements  since  1966. 
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The  literature  review  of  those  articles  written  before  1966  has 
been  completed,  and  the  findings  are  reviewed  herein.  Preliminary  to 
the  actual  review  is  a  description  of  the  situation  which  is  under  study. 
The  appropriate  conservation  laws  are  developed  but  cannot  be  solved 
for  the  properties  of  the  combustion  system  since  the  analysis  is  mathe¬ 
matically  intractable.  Thus,  it  must  be  replaced  by  an  approximate 
method  which  relies  on  the  theory  of  superpositioning  as  utilized  in 
the  Northern  Research  analysis. 

The  remainder  of  the  report  consists  of  a  review  of  the  analytical 
and  empirical  techniques  available  for  prediction  of  wall  temperature 
distributions.  As  noted  previously,  each  of  the  articles  reviewed  was 
written  before  1966.  The  three  basic  modes  of  heat  exchange,  convection, 
conduction,  and  radiation  are  covered  in  the  forms  found  in  the  combustion 
chamber.  The  final  section  of  the  report  consists  of  a  sumnary  of  the 
analyses  presented  and  some  conclusions  as  to  the  direction  of  advance¬ 
ment  since  1966. 

Before  embarking  on  the  development  of  a  suitable  model,  the  connec¬ 
tion  between  this  portion  of  the  gas  turbine  program  and  that  discussed 
in  Part  I  (Hammond  and  Mellor,  1971)  and  Part  III  (Anderson  and  Mellor, 
1971)  of  the  Second  Annual  Report  should  be  made  clear.  The  analysis 
of  Hammond  and  Mellor  (1970a,  1970b)  provides  the  necessary  mean  gas 
temperature  and  composition  along  the  centerline  of  the  combustor  for 
any  heat  transfer  analysis.  Another  connection  between  the  combustor 
modelling  project  and  the  heat  transfer  analysis  is  the  NREC  airflow 
and  heat  transfer  analysis  (Anon.,  1968),  which  has  been  incorporated 
in  a  computer  program  written  for  NASA  Lewis  Research  Center.  This  pro¬ 
gram  provides  the  air  flow  distribution  necessary  for  the  determination 
of  the  mean  temperature  and  gas  composition  from  the  combustor  modelling 
program  and  as  noted,  its  heat  transfer  program  provides  a  wall  tempera¬ 
ture  distribution. 

Inmediately  the  question  arises  as  to  why  the  NASA  Lewis  program 
cannot  be  used  for  the  heat  transfer  analysis.  The  reason  is  that  to 
determine  the  wall  temperature  distribution  with  this  program,  the  oper¬ 
ator  must  choose  which  modes  of  heat  transfer  are  most  important  in  the 
combustor  chamber  he  is  designing.  The  only  means  available  at  the 
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present  time  for  determining  whether  the  correct  choice  has  been  made 
Is  to  compare  the  analytical  results  with  experimental  data  for  the 
combustor.  Therefore,  the  program  Is  of  limited  use  since  empirical 
determinations  of  heat  transfer  are  still  required.  Clearly  it  is  de¬ 
sirable  to  develop  a  purely  analytical  technique  for  heat  transfer 
estimates,  and  the  present  study  Is  concerned  with  determining  if  this 
possibility  exists  In  view  of  the  literature  published  since  1966. 

In  order  to  ascertain  the  validity  of  any  advancements  made  in 
the  analytical  determination,  the  experimental  results  for  the  wall 
temperature  distribution,  to  be  obtained  in  the  facility  described  in 
the  final  volume  of  this  report  (Anderson  and  Me 11 or,  1971)  will  be 
necessary. 

In  summary,  the  purpose  here  is  to  report  on  the  first  phase  of 
a  literature  review  on  heat  transfer  in  gas  turbine  combustion  chambers. 
The  present  report  Is  limited  to  the  pre-1966  period  and  thus  allows 
appreciation  of  the  state  of  the  art  at  the  time  of  writing  of  the  NASA 
Lewis  program,  as  well  as  advances  available  in  the  literature  since 
1966.  A  later  report  will  be  concerned  with  post-1966  literature  and 
the  feasibility  of  Improving  existing  heat  transfer  analysis  of  gas  tur¬ 
bine  combustors. 


CHAPTER  II 


THE  MODEL 


A.  Basic  Equations 

The  purpose  of  this  section  is  to  establish  the  analytical  methods 
available  for  modelling  heat  transfer  in  the  combustion  chamber.  The 
physical  phenomena  of  interest  occurring  in  the  chamber  are  defined  by 
the  turbulent  flow  conservation  equations.  The  presence  of  recircula¬ 
tion  and  the  rotational  dependence  of  the  flow  will  be  shown  to  make 
this  approach  useless.  The  final  section  reviews  the  approximations 
that  have  been  made  to  analytically  account  for  the  heat  transfer  to 
the  combustor  wall. 

Figure  1  illustrates  a  typical  combustor  of  turbular,  tubo-annular , 
or  annular  design.  For  each  region  of  the  chamber  presented,  the  primary 
zone,  the  seconary  zone,  the  annulus  flow,  and  the  chamber  wall,  the 
laws  of  conservation  of  mass,  momentum,  species,  and  energy  are  appli¬ 
cable.  Each  control  volume  shown  for  the  four  zones  (C.V.l,  C.V.2,  C.V.3, 
and  C.V.4)  requires  a  somewhat  different  form  of  the  conservation  equa¬ 
tions.  The  equations  are  simply  presented  here  because  very  adequate 
derivations  exist  in  numerous  texts  (Schlichting,  1960,  Bird,  Stewart, 
and  Lightfoot,  1960). 

Conservation  of  Mass 
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Conservation  of  Species 
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Conservation  of  Energy 
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A  nomenclature  list  is  provided  at  the  beginning  of  the  report.  The 
| above  equations  are  for  steady,  incompressible,  turbulent  flw. 

The  continuity  equation  for  steady  inconpressihle  flow  is  self- 
; explanatory.  In  the  momentum  equation,  the  momentum  diffusion  term, 
i a/3xi  [(u  +  A  )  3ulj/3x. 3  accounts  for  the  laminar  viscosity  as  well  as 
ithe  Reynolds  stress  term  3/3Xf[Ax  3^/ax.].  The  rate  of  gewation  or 
destruction  of  species  m  is  designated  by  the  last  term  in  Eqn.  (2-3). 
*For  the  energy  equation,  the  term  on  the  left  hand  side  represents  the 
convective  transfer  of  thermal  energy.  The  first  term  on  the  right 
hand  side  represents  the  laminar  thermal  diffusion  plus  the  eddy  or 
apparent  dlffusivity  due  to  turbulence  (Schlichting,  1960).  The  radia¬ 
tion  exchange  and  the  heat  generation  per  unit  volume  are  represented 
by  the  next  tow  terms  on  the  right  side.  The  thermal  expansion  pressure 
and  the  viscous  dissipation  terms  are  the  final  terms  in  this  equation. 
The  dissipation  term,  $,  represents  the  turbulent  viscous  heating  for 
which  the  coefficient  (u  +  A  )  represents  the  laminar  and  turbulent 
[(eddy)  viscosities. 

The  bars  over  the  various  components  indicate  the  time  average 
property  of  the  particular  term.  Each  term  consists  of  a  time  average 
component  and  a  component  due  to  the  velocity  of  fluctuation.  To  cite 
an  example,  the  velocity  u^  is  given  by; 


v/here  u^  is  the  time  average  velocity  and  u-  is  the  fluctuation  velocity. 
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This  term  is  incorporated  in  the  dissipation  function  $  and,  for 
momentum,  in  the  3tTj/3x-  term  associated  with  viscosity  (Schlichti ng, 
1960). 

The  four  control  volumes  depicted  In  Fig.  1  are,  in  general,  de¬ 
scribed  completely  by  Eqns.  (2-1)  -  (2-4)  plus  an  equation  of  state. 
However,  for  all  but  one  of  these  control  volumes,  the  primary  zone, 
it  is  possible  to  eliminate  terms  or  one  of  the  equations  and  still 
determine  all  unknowns.  The  process  of  eliminating  terms  will  provide 
useful  physical  insight  into  the  phenomena  occurring  in  the  zones  where 
each  control  volume  is  situated. 

For  the  primary  zone  control  volume,  the  process  of  combustion  re¬ 
quires  all  four  of  the  conservation  laws  and  an  equation  of  state  to 
determine  the  unknowns  since  combustion,  and  therefore  a  change  in 
species,  requires  the  use  of  the  species  equations  to  determine  mass 
fractions  of  those  species  present.  Therefore,  for  the  primary  zone 
all  of  the  equations  in  the  form  stated  are  necessary. 

As  a  first  approximation,  in  the  secondary  zone  the  process  of  com¬ 
bustion  can  be  considered  to  have  either  gone  to  completion  or  termin¬ 
ated.  The  species  concentrations  are  unchanged  by  chemical  reaction 
but  still  vary  due  to  the  secondary  air  entrainment.  Although  the 
species  equation  must  be  used,  it  reduces  to 
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The  cessation  of  the  chemical  processes  also  allows  the  elimination  of 
the  heat  generation  term,  q^",  in  the  energy  conservation  equation, 
Eqn.  (2-4). 

The  final  gas  phase  control  volume,  the  annulus  volume,  C.V.3,  is 
taken  to  have  a  constant  mass  fraction  of  oxygen  and  nitrogen  since 
combustion  does  not  occur  here.  The  species  equation  and  heat  genera¬ 
tion  term  may  thus  be  neglected.  Stemming  from  the  low  temperatures 
in  the  region  and  the  knowledge  that  nitrogen  and  oxygen  do  not  part¬ 
icipate  in  radiation  transfer  (McAdams,  1954),  the  radiation  transfer 
rate  need  not  be  considered  in  Eqn.  (2-4).  The  existence  of  only  two 
species  has  truncated  the  system  of  equations  to  four,  Eqns.  (2-1), 
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(2-2),  (2-4)  and  an  equation  of  state. 

The  final  region  of  interest  is  the  solid  phase  combustor  wall 
or  liner.  Assuming  the  control  volune  in  the  wall  is  of  infinitesmal 
axial  length  but  of  finite  thickness  (equal  to  the  wall  thickness) 
will  simplify  the  solution  to  the  equations.  The  wall  is  assumed  to 
be  stationary,  of  constant  mass  and  chemical  composition,  and  has  no 
internally  stored  stresses  or  internal  heat  sources.  If  in  addition 
body  forces  are  neglected,  Eqns.  (2-1)  -  (2-3)  are  eliminated  from  the 
problem,  and  the  energy  equation  is  given  by 
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For  most  materials  the  thermal  expansion  term  is  small  (Guy,  1959)  and 
the  material  is  opaque  to  radiation.  Eqn.  (2-6)  may  be  further  simplified 
by  eliminating  the  turbulent  conduction  term,  which  yields  the  final 
equation: 


Under  the  assumption  that  the  temperature  is  uniform  normal  to  the  wall, 
(in  the  x?  -  direction)  the  equation  may  be  easily  expanded. 


^  ^  +  isj  ^3  ■  0  (2-8) 

Note  that  the  turbulent  temperature  notation  has  been  removed  because 
the  concept  of  a  turbulent  temperature  in  a  stationary  solid  has  no 
meaning.  One  further  assumption  which  has  been  made  by  Winter  (1955), 
tefebvre  and  Herbert  (1960),  and  Tipler  (1955)  and  followed  by  NREC 
(Anon.,  1968),  is  that  conduction  in  the  tangential  direction,  x3,  is 
neglected.  Thus,  Eqn.  (2-8)  yields: 
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The  effects  of  radiation  heating  or  cooling  and  convection  appear  to 
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be  neglected  In  this  equation.  However,  they  are  found  to  be  Included 
In  the  boundary  conditions  for  the  chamber  wall  control  voli**. 

In  conclusion,  the  equations  necessary  for  a  solution  of  the  flow 
and  temperature  fields  become  Increasingly  more  simplified  as  the 
phenomena  occurring  In  the  region  decrease  In  mmfcer.  The  primary  zone 
requires  the  solution  of  a  very  complicated  set  of  equations  {Eqn.  (2-1)- 
(2-4))  Involving  reaction  kinetics.  Cessation  of  combustion  simplifies 
the  system  of  equations  for  the  secondary  zone  by  eliminating  the  chemical 
kinetics  and  the  term  representing  an  Internal  heat  source  In  the  energy 
equation  (Eqn.  (2-4)).  Annulus  flow  does  not  require  the  above  terras 
and  laws  the  additional  exclusion  of  radiation  participation.  Finally, 
the  sulirf  contr.  :  volune  requires  only  Eqn.  (2-9)  and  two  boundary 

conditions.  The  complexity  of  the  flow  makes  the  solution  for  the  gas 
phase  prohibitive,  and  only  the  solid  phase  control  volume  will  be  con¬ 
sidered  beyond  a  few  concluding  remarks. 

An  experimental  Investigation  by  Hlett  and  Powell  (1962)  indicated 
the  recirculation  zone  found  In  a  typical  combustor  (Fig.  2).  Coupling 
fluid  motion  such  as  this  with  the  Induced  rotational  motion  due  to  the 
swlrlers  which  are  present  In  most  combustors  makes  the  solution  of  Eqns. 
(2-1)  -  (2-4)  impossible.  For  this  reason,  other  than  a  theoretical 
heat  release  rate,  the  conditions  In  the  gas  phase  have  been  eleminated 
from  the  analyses  of  Winter  (1955),  Lefebvre  and  Herbert  (1960),  Tipler 
(1955),  and  NR EC  (Anon.,  1968). 

Therefore,  the  remainder  of  this  chapter  will  review  the  approxi¬ 
mations  made  In  the  literature  to  replace  the  solution  for  the  gas  phase 
properties.  The  mathematically  expedient  process  of  super  positioning 
will  be  shown  to  approximate  the  solution  for  the  wall  temperature  dis¬ 
tribution  and  thus  avoid  the  nonlinear  Integro-differential  equations 
necessary  for  a  mathematically  exact  solution.  Utilizing  superposltlon- 
Ing  techniques  will  also  separate  the  Individual  components  of  heat 
transfer  into  six  terms,  which  may  then  be  discussed  separately. 

B.  Interaction  of  Reflation,  Convection, 
and  Conduction 

Determination  of  the  energy  f luxes  to  the  wall  requires  an  under- 
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standing  of  radiation  transport  in  a  participating  medium  and  its  in¬ 
teraction  with  convective  and  conductive  heat  transfer.  A  prerequisite 
is  thus  a  development  of  the  radiative  theory  of  transport.  Viskanta 
(1966)  developed  the  equation  of  transport  on  a  macroscopic  basis  and 
this  work  shall  be  reviewed.  Utilizing  a  similar  development,  Einstein 
(1963)  developed  the  equations  necessary  for  determining  the  temperature 
field  in  a  flow  of  air  through  a  passage  including  convective  and  con¬ 
ductive  terms  as  well  as  radiation  participation.  The  result  of 
Einstein's  work  is  the  substantiation  of  the  theory  of  superpositioning. 
Superpositioning  of  energy  transfer  modes  implies  that  radiation,  con¬ 
vection,  and  conduction  are  each  independent  of  the  other  two  contribu¬ 
tions.  The  final  topic  of  this  section  will  be  the  development  of  the 
energy  exchange  equation  for  the  wall  using  the  theory  of  superpositioning. 
The  exchange  equation  will  be  the  basis  for  discussing  the  literature 
written  before  1966. 

Viskanta  (1966)  has  provided  a  derivation  of  the  radiation  trans¬ 
port  theory  based  on  the  physical  interaction  of  radiation  witn  a  volume 
of  gas.  Two  directions  are  possible  for  the  development  of  the  trans¬ 
port  theory:  the  microscopic  and  macroscopic  theories.  Microscopic 
theory  concerns  the  interaction  of  electromagnetic  waves  with  the  parti¬ 
cular  species  present  in  the  medium.  The  macroscopic  theory  accounts 
for  the  emission,  absorption,  or  scattering  of  radiation  upon  inter¬ 
action  with  a  volume  of  gas.  Emphasis  must  be  placed  on  the  fact  that 
radiation  is  attenuated  by  matter  and  not  by  a  volume;  the  study  on  a 
macroscopic  basis  is  thus  the  overall  interaction  of  radiation  with 
the  matter  present.  In  the  following,  the  macroscopic  case  has  been 
considered.  The  assumptions  Viskanta  used  are  reviewed  and  the  transport 
equation  stated  with  a  minimum  of  proof. 

Several  assumptions  are  necessary  to  formulate  a  problem  which  is 
mathematically  tractable.  The  theory  of  radiation  transport  ignores  the 
electromagnetic  wave  phenomena  of  i nterference ,  diffraction,  and  coherence, 
and  is  thus  a  limiting  case  of  the  more  general  electromagnetic  theory 
of  radiation  transport.  Since  in  order  to  determine  the  variation  in 
temperature  in  the  medium  it  is  necessary  to  define  a  temperature ,  it 
is  assumed  that  local  thermodynamic  equilibrium  exists.  Thus,  non- 
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equilibrium  radiation  such  as  chemiluminescence  is  neglected.  The 
third  assumption  made  is  that  the  medlun  is  continuous  and  homogeneous 
and  is  capable  of  absorbing  and  emitting  radiation  in  a  coherent,  iso¬ 
tropic  manner.  That  is,  radiation  scattered  by  the  matter  in  the 
medium  does  not  have  a  preferred  direction  and  remains  at  the  frequency 
of  Incidence. 

Consider  now  a  cylinder  of  participating  media  which  is  subjected 

to  a  pencil  of  radiation  entering  through  one  end  (Fig.  3).  The  ray  is 

travelling  in  a  direction  s  confined  within  a  solid  angle  dQ.  The 

monochromatic  intensity,  Iv,  entering  the  control  volume  is  attenuated 

such  that  the  intensity  at  the  end  of  the  control  volume  is  I  +  dl  . 

v  v 

The  change  in  intensity  is  determined  by  the  radiation  transport  equation. 

The  radiation  attenuation  by  absorption,  scattering,  and  emission 
may  be  expressed  mathematically.  The  spectral  intensity  of  radiation 
at  the  end  of  the  cylinder  may  be  expressed  as: 

I  dA  dQ  dv  dt 
v 

Inside  the  volume,  energy  may  be  lost  from  the  pencil  due  to  absorption 
and  scattering,  which  may  be  represented  as  follows 

(a  +  k  )  I  ds  dA  dv  dt  dQ 
v  v'  v 

Note  the  inclusion  of  the  elemental  distance  ds  traveled  along  s.  The 
two  coefficients  are  the  coefficients  of  monochromatic  scattering 
and  absorption  The  matter  within  the  volume  emits  radiation  in  a 
frequency  range  from  v  to  v  +  dv  in  the  direction  of  s,  and  this  energy 
may  be  represented  by: 

dV  dQ  dv  dt 

The  radiant  intensity  is  increased  by  radiation  scattered  into  the 
volume  from  all  possible  directions  which  is  represented  by: 
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ev  ds  CO/47T)  /n,=47r  pv(s*  -  s)  Iv  (s')  da']  dA  da  dv  dt 


The  function  p  (s'  -  s)  is  related  to  the  probability  of  a  photon  en- 
tering  the  volume  traveling  along  s'  and  being  monochroma tically  scat¬ 
tered  so  that  its  final  direction  is  s.  The  scattering  function  has 
been  normalized  to  unity  since  the  probability  over  all  directions  must 
be  unity: 


i/4ir/Ql=4ii  pv(s‘  -  s)  da'  =  1 


(2-10) 


The  attenuation  due  to  absorption,  emission,  and  scattering  may 
be  formulated  into  an  equation  representing  the  conservation  of  radiant 
intensity  over  the  elemental  length  ds  in  a  frequency  range  v  to  v  +  dv, 
in  a  time  dt  and  area  dA: 


(Iv  +  dlv  -  Iv)  dA  da  dt  dv  =  dly  dA  da  dt  dv 


(2-11) 


The  attenuation  dlv  dA  da  dt  dv  can  be  represented  by  those  losses  and 
gains  mentioned  above. 


dlv  dA  da  dv  dt  =  nv  dV  da  dv  dt  -  (a^  +  <y)  Iv  ds  dA  da 
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Iv  (s')  da']  dA  da  ds  dv  dt 


(2-12) 


Since  dV  equals  dA  ds,  this  equation  may  be  divided  by  dA  ds  da  dv  dt. 
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Noting  that  the  distance  ds  may  be  represented  by  c  dt: 
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where 
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In  actuality,  the  derivative  of  I  with  respect  to  time  may  be  replaced 
by  the  substantial  derivative: 
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where,  by  definition: 
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Therefore , 
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For  most  engineering  applications  the  time  derivative  is  negligible 
because  of  the  1/c  factor.  Thus, 


n-(a  +k)I  +  a  H 
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(2-18) 


which  Is  known  as  the  quasi-steady  radiation  transport  equation. 

The  conservation  of  radiant  energy  may  be  formulated  using  the 
above  equation  by  integrating  over  all  solid  angles  (£2=4tt): 
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where  6  is  the  monochromatic  extinction  coefficient  and  is  defined  as 
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energy  equation  is  the  fact  that  scattering  has  been  eliminated  from 
the  equation.  This  implies  that  scattered  energy  is  not  stored  by  the 
medium  which  is  physically  true. 

The  conservation  of  radiant  energy  equation  may  be  related  to  the 
previously  formulated  energy  equation  (2-4)  by  replacing  the  3q"*/3x^ 
term  by  a^/9x^  where  7\j  is  the  component  of  radiative  flux  in  the  i 
direction  integrated  over  all  frequencies.  Thus, 


Integrating  over  all  frequencies  and  noting  that  the  emission  for  a 
medium  in  local  thermodynamic  equilibrium  is  given  by: 


VT>  =  \  hv  <T>  '2-?7> 

where 

3 

2h  v 

'»v (T)  ■  (2-28) 

c0  (e  H  -  1) 


yields 


C  %  (T)  dv  -  C  \  Jbv  <T>  *■ 


(2-29) 


Utilizing  these  equations  for  the  variation  in  radiant  flux  along  the 
x..  direction: 

77  aT 


aT 


p  cp  ui  axi  =  axi  ^Ki  +  Ag  cp  ax  J  "  axi  ^  £  Kv  !bv  ^  * 


dv  -  /  Kv 


dn  dv]  +  q" 1  +  bT  |£- 


+  (u  +  At)  4) 


(2-30) 
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Aside  from  the  very  difficult  problem  of  determining  the  necessary 
coefficients  for  the  above  equation,  a  mathematical  solution  in  closed 
form  would  be  extremly  difficult  to  obtain  except  for  the  simpliest 
cases. 

Einstein  (1963)  noted  that  the  effects  of  radiation  and  conduction 
could  be  accounted  for  by  superimposing  the  energy  transport  due  to 
each  of  the  modes  in  the  absence  of  the  other.  The  system  considered 
is  shown  in  Fig.  4.  The  effects  of  absorption  of  the  gas  in  such  a  con¬ 
vective,  conductive  system  was  studied  utilizing  an  equation  very  simi¬ 
lar  to  Eqn.  (2-30)  but  for  laminar  flow: 


4k  aSB  T4  »)  +  D  «  cp  fl&  /;4()  -  K  /f4o 


=  k  Iff  aSB  T4  (r )  f (?'-  ftn)  dt 


+  k  //  a$B  Ts4  (r)  g(r  -  ftQ)  dA  (2-31)* 

i  ; 

The  definitions  of  the  individual  terms  as  presented  by  Einstein  are 

i  f 

given  in  Table  1.  Einstein  used  Eqn.  (2-31)  to  determine  the  tempera- 

;  j  | 

ture  distribution  in  the  channel  for  a  constant  wall  temperature.  In 

j  ♦  ,f 

order  to  solve  the  integro-differential  equation,  it  was  necessary  to 

j  ’  i 

replace  the  differentials  with  approximate  finite  differences  in  tem¬ 
pi 

perature  between  discrete  zones  in  the  flow.  A  grid  of  one  hundred 

units  was  used.  In  this  manner  he  was  able  to  determine  the  heat 

l  I 

transfer  between  the  two  plates  and  the  temperature  profile  between 
the  plates. 

Upon  analyzing  the  effects  of  radiation  and  conduction  in  the 
absence  of  convection,  the  total  heat  transferred  was  slightly  greater 
than  the  heat  transferred  by  conduction  in  the  absence  of  radiation 
summed  with  the  heat  transferred  by  radiation  in  the  absence  of  conduc¬ 
tion.  This  implies  that  the  total  heat  transfer  may  be  approximated 


★  .  . 

Note:  The  nomeclature  for  Eqn.  (2-31)  is  defined  in  Table  1. 
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Table  1. 

Terms  In  the  Einstein  (1963)  Energy  Equation 


Term 


4k  oSB  T-(R0) 


PUC 

M  p  3  X 

3y  r  o 

k  ///  aSB  T4  (f)  f  (r  -  f$Q)  dt 


k  //  a$B  Ts4  (r)  g  (r  -  $Q)  dA 


_ Description _ 

Radiant  energy  emitted  per  unit 
voline  at  r  =  RQ 

Enthalpy  increase  per  unit  volume 
of  the  flowing  gas  at  ?  =  RQ 

Net  conduction  heat  transfer 
into  the  unit  volume 

Radiation  absorbed  per  unit 
volume  at  R0  from  emission  given 
off  by  surrounding  gas 

Radiation  absorbed  per  unit  volume 
at  R0  from  emission  of  flat  plates 
and  end  surfaces 


-Zi  nstei  n '  s  Nomenc  lature 

k  =  radiation  absorbtion  coef.  of  gas 

ft  =  position  vector  of  variable  point  in  channel  or  surface 

r  =  position  vector  of  variable  point  in  channel  or  surface 

dx  s  infinitesimal  volume  element  in  gas 

f(s)  =  gas-1 ine -source  to  gas  radiation  exchange  factor 

g(s)  =  surface-line  source  to  gas  radiation  exchange  factor 

Subscripts 
s  =  surfacl 

o  =  integrlted  mean  conditions  at  channel  outlet 
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by  the  summation  of  the  individual  components  acting  independently  of 
each  other: 


qRC  *  qR  ♦  qc  (2-32) 

The  accuracy  of  this  assumption  is  exhibited  in  Fig.  5;  qRC/(qR  +  Qq) 
is  plotted  against  the  optical  thickness  rQ  for  various  values  of 
NrR  =  K/Dp  aSB  T  ,  which  is  a  measure  of  the  importance  of  conduction 
coriipared  to  radiation  in  the  system. 

Einstein's  results  indicate  that  the  approximation  is  most  accurate 
for  the  two  limits  of  an  optically  thick  or  thin  medium  (the  optical 
thickness  is  a  measure  of  the  mean  free  path  of  a  photon  in  the  system). 
As  would  be  expected,  for  heat  transfer  when  conduction  (TLR  =  2.08) 

or  radiation  (NCR  =  0.0208)  is  dominant,  the  approximation  has  a 
maximum  error  of  4%.  When  both  modes  influence  the  heat  transfer  (N^R 
=  0.104  or  fJCR  0.208)  the  error  is  a  maximum  at  9%.  Viskanta  and 
Grosh  (Cess,  1964)  utilizing  a  similar  physical  system  found  that  the 
maximum  error  was  11%.  Again,  the  discrepancy  was  greatest  for  systems 
which  had  significant  radiation  and  convection. 

An  extension  of  the  above  example  of  superpositioning  has  been 
utilized  by  Winter  (1955),  lefebvre  and  Herbert  (1960),  and  Tipler  (1955) 
Rather  than  considering  the  system  of  equations  (2-1)  -  (2-4)  for  each 
of  the  gas  control  volumes  to  determine  the  boundary  conditions  for  Eqn. 
(2-9),  the  solid  control  volume  equation,  the  process  of  heat  transfer 
has  been  approximated  as  the  sunmation  of  the  individual  components  of 
heat  transfer  acting  as  if  the  other  components  were  not  present. 

Fig.  6  illustrates  the  analysis  of  Tipler  (1955)  using  the  concept 
of  superpositioning.  The  radiant  heat  flow  from  the  flame  to  the  com¬ 
bustion  chamber  wall  is  represented  by  R^.*  The  radiant  exchange  bet¬ 
ween  the  walls  at  high  temperature  (in  the  primary  zone,  for  example) 
to  the  walls  at  lower  temperatures  (in  the  secondary  zone)  is  designated 
by  R^. j  The  chamber  walls  will  also  increase  in  temperature  due  to  the 

i  f  I 


•fc  : 

The  notation  used  herein  is  attributed  to  Winter  (1955). 
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MODES  OF  HEAT  TRANSFER 


R]  +  C]  +  R3  =  R2  +  C2  +  K 


FIGURE  G  TIPLER'S  MODEL  OF  HEAT  TRANSFER  COMPONENTS 
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convective  heating  of  the  hot  combustion  gases  represented  by  C-j. 

Assuming  for  the  moment  that  the  wall  is  insulated,  the  heat  transfer 
to  a  wall  element  may  be  calculated  in  the  following  manner: 

K  +  R1  +  C]  =  R3  (2-33) 

The  sign  of  K,  the  conduction  component,  is  dependent  on  the  gradient 
of  the  temperature  as  defined  by  Fourier’s  law  of  heat  conduction 
(McAdams,  1954): 

K  =  Ki|L-  (2-34) 

In  the  instance  of  secondary  air  flow  through  the  annulus,  heat 
will  be  lost  from  the  wall  by  convection  and  radiation  (Tipler,  1955) 
since  the  anr.ulus  gases  are  at  a  lower  temperature  than  the  combustion 
products  in  the  chamber.  The  convective  cooling  component  of  the 
annulus  gases  is  specified  by  C^.  Three  other  convective  cooling  methods 
could  be  present  in  a  typical  combustor:  penetration  jets,  film  cooling 
slots,  and  transpiration  cooling.  At  present  these  deviations  from 
the  flat -plate -type  convective  cooling  will  be  neglected.  Radiation 
exchange  R^  is  also  present  in  the  annulus  in  the  form  of  radiation 
transfer  between  the  chamber  wall  and  the  chamber  casing.  The  annulus 
terms  may  be  considered  losses  from  the  wall  element: 


R-|  +  C-j  +  K  =  R3  +  ^2  +  ^2 


(2-35) 


The  remainder  of  this  report  deals  with  the  investigation  of  each 
of  the  terms  in  Eqn.  (2-35),  and  the  permutations  found  in  typical  com¬ 
bustors.  In  the  next  chapter,  radiation  in  the  chamber  and  annulus  is 


considered,  and  the  appropriate  methods  presented  in  the  literature  of 
interest  are  reviewed,  for  the  most  part,  in  an  order  of  increasing 
approximation  or  dependence  on  empirical  results.  In  Chapter  IV  the 


theories  available  for  predicting  the  magnitude  of  each  of  the  convective 

•> 

terms  are  presented.  Throughout:  the  report,  discussion  of  empirical 

i  • 

results  has  been  kept  to  a  minimum  and  the  emphasis  is  placed  on 
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theoretical  developments,  since  an  analytical  model  is  desired. 


CHAPTER  III 


RADIATION  TRANSFER  IN  COMBUSTION  CHAMBERS 


A.  Introduction 


Radiation  transfer  In  furnaces  has  been  studied  on  an  empirical 
and  limited  theoretical  basis  for  many  years  (Lefebvre,  1965).  How¬ 
ever,  only  since  1960  have  the  theories  of  absorption  and  emission  by 
a  medium  been  applied  to  engineering  rather  than  astrophysics  problems 
(Cess,  1964).  The  result  of  this  theory,  the  equations  of  radiation 
transport  (Eqn.  (2-18))  and  the  conservation  of  radiant  energy  (Eqn. 
(2-25)),  are  unfortunatly  restricted  by  their  mathematical  coirplexity. 
This  leads  directly  to  the  discussion  of  two  approximate  methods,  the 
one  dimensional  and  two  dimensional  formulae  for  flame  radiation.  Also 
necessary  for  the  determination  of  flame  radiation  is  a  method  for 
calculating  the  emisslvitles  and  absorptivities  of  the  flame.  The 
correlations  available  for  these  calculations,  for  both  luminous  and 
nonluminous  flames,  will  be  covered  in  the  second  section.  The  final 
two  topics  covered  in  this  chapter  are  the  evaluation  of  and  R^, 

the  exchange  between  the  wall  of  the  chamber  and  the  casing  and  the 
exchange  between  different  surface  areas  on  the  inside  of  the  chamber, 
respectively . 


B.  Flame  Radiation 


Of  the  work  related  to  the  determination  of  heat  transfer  in  the 
combustion  chamber,  a  great  deal  has  been  devoted  to  the  study  of  radiant 
exchange  between  the  flame  and  surroundings  (Lefebvre  and  Herbert,  1960). 
Unfortunatly,  a  large  portion  of  this  work  concerns  the  development  of 
experimental  techniques  for  investigating  this  phenomenon.  The  two 

methods  of  calculating  the  magnitude  of  R-,  which  have  been  presented  in 

i  f  1 

the  literature  depend  substantially  on  empirical  results  for  accuracy. 
Both  of  the  formulations  and^  their  restrictions  are  discussed  in  detail 
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below.  The  final  topic  of  Interest  In  flame  radiation  Is  an  under¬ 
standing  of  the  empirical  results  and  correlations  of  these  results 
used  to  determine  flame  emisslvities  for  both  nonluminous  and  luminous 
flames. 

The  one  dimensional  exchange  of  radiation  between  the  flame  and 
the  wall  at  the  same  axial  position  as  the  flame  was  first  applied  to 
aircraft  combustion  chambers  by  Winter  (1955).  Emphasis  must  be  placed 
on  the  fact  that  the  analysis  is  applicable  to  two  dimensional  flame 
radiation  (that  is,  radiation  exchange  between  the  flame  element  and 
the  wall  at  all  axial  positions)  only  under  the  assumption  of  a  constant 
wall  temperature  (McAdams,  1954).  For  the  case  considered  here,  the 
exchange  is  limited  to  elements  of  the  flame  and  isothermal  elements 
of  the  wall  at  the  same  axial  position.  If  sufficiently  small  wall 
elements  are  considered  the  area  is  indeed  isothermal. 

Utilizing  the  formulation  for  the  thermal  radiant  emission  from 
a  black  body: 


(3-1) 


the  energy  arriving  at  the  wall  due  to  radiation  is  proportional  to 
the  difference  between  the  energy  emitted  from  the  flame  and  the  energy 
radiated  from  the  wall  which  is  absorbed  by  the  flame. 

R1  *  ^ef  °SB  Tf  "  af  CSB  Tw  ^ 

The  coefficients  and  represent  the  emissivity  and  absorptivity 
of  the  flame,  respectively.  Eqn.  (3-2)  is  derived  for  the  specific 
case  of  a  chamber  wall  which  is  black,  c,  =  1.0.  To  account  for  walls 
which  are  not  black  and  reflect  radiation,  an  additional  factor  o' 
has  been  incorporated  into  the  equation  (McAdams,  1954).  The  factor 
accounts  for  radiant  energy  which  after  leaving  the  flame  has  been 
reflected  back  through  the  flame  and  is  absorbed  or  rereflected  by  the 
wall  on  the  opposite  side  of  the  flame.  Empirical  results  (McAdams,  1954, 
Winter,  1955)  have  indicated  that: 
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satisfactorily  accounts  for  the  multiple  reflections  encountered  in 
the  chamber. 

Lefebvre  and  Herbert  {I960)  have  found  that  over  a  wide  range  of 
optical  thicknesses  the  absorptivity  and  emissivity  of  the  flame  are 
related  by  the  following 


Tf  1.5 
(/> 


w 


(3-4) 


In  the  limit  of  flame  and  wall  temperature  equality,  af  =  as  would 
be  expected  from  Kirchoff's  law  (McAdams,  1954).  Eqn.  (3-3)  and  (3-4) 
may  be  substituted  into  the  equation  for  the  flame  radiation  to  give: 


1  +  c 

'  (  1  )  ef 


SB  *f 


tJ.5  (t  2.5  _  Tw2.5) 


(3-5) 


which  is  explicitly  valid  for  the  one  dimensional  case  or  two  dimensional 
case  under  the  restrictions  outlined  previously. 

Two  values  are  still  necessary  for  the  calculation  of  R-j.  In  the 
above  formula,  radiant  energy  does  not  appear  to  be  dependent  on  the 
geometry  of  the  system.  However  the  geometric  factor  is  found  in  the 
determination  of  flame  emissivities,  which  will  be  discussed  in  detail 
after  the  second  method  of  calculation  is  reviewed.  A  second  restriction 
to  immediate  calculation  is  the  determination  of  the  flame  temperature. 

A  mean  temperature  for  the  flame  is  usually  determined  from  the  temper- 
aturerisedue  to  combustion  over  the  inlet  temperature,  TJN: 

Tf  =  Tin  +  AT  (3-6) 


The  value  of  AT  may  be  determined  from  standard  temperature  rise  curves 
based  on  the  local  fuel  to  air  ratio  (Lefebvre  and  Herbert,  1960). 

Neglecting  the  radiation  from  a  flame  element  to  a  section  of  the 
wall  other  than  the  section  at  the  sane  axial  location  has  obvious 
physical  disadvantages.  For  this  reason,  Hottel  (McAdams,  1954) 
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extended  the  one  dimensional  analysis  to  account  for  exchange  between 
the  flame  and  all  sections  of  the  wall  and  removed  the  restriction  of 
constant  temperature  along  the  wall.  Hottel  began  by  specifying  the 
analysis  necessary  for  determining  the  shape  of  the  gas  volume. 

If  a  hemisphere  of  gas  is  considered  to  be  radiating  at  a  rate 

4 

a$B  £q  Tq  to  a  surface  element  at  the  base  of  the  hemi sphere ,  the 
emission  from  an  element  at  rad1usrg  to  the  element  per  unit  area,  per 
unit  of  shell  thickness  is  {see  Fig.  7): 


d(°SB  CG  O 


*  aSB  TG 


4  deG 

*7 


(3-7) 


The  Irradiation  incident  on  the  surface  per  unit  thickness  of  the  shell 
dr$,  is  the  radiant  emission  divided  by  a  factor  of  it: 


dl  -  [- 


'SB  'G 


dcr 

aiy1  drs 


(3-8) 


t 


Utilizing  the  differential  exchange  between  the  volume  and  the  elemental 
surface  area, 


dR-j  *  dl  dA  cos  dfl 


SB 


V  der 

■ - ]  dA  dfl  cos  $  gy  dr$ 


(3-9) 


Integration  over  all  of  the  volume  and  for  the  base  area  yields: 

R1  *  'a  'o"“  Wos  *  dn  drs  "  (3'10) 

The  Integration  of  the  above  is  a  tedious  procedure  due  to  the  required 
determination  of  cos  <j>  dn  for  the  geometry  considered.  A  similar  and 
simpler  procedure  was  employed  to  determine  the  mean  beam  length  I 
such  that  the  radiant  energy  emitted  over  the  entire  volume  is  equal 
to  the  average  rate  of  emission  times  the  mean  beam  length.  The  ^alue 
of  I  was  found  to  depend  on  the  partial  pressure  of  the  participating 
gas,  Pg,  times  a  characteristic  length  X.  This  dependence  is  almost  a 
constant  function  of  the  value  of  X,  which  has  been  determined  for 
certain  geometries  of  interest  (McAdams,  1954). 

Hottel  has  shown,  utilizing  Eqn.  (3-7),  that  the  emission  per  unit 
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volume  of  gas  Is: 


(a-TTr)/  - 

d  PG  Pq£=0 


(3-11) 


Tha  value  of  the  differential  of  the  gas  emissivity  is  evaluated  as  if 
no  radiation  is  absorbed  by  the  intervening  gases.  The  partial  pressure 
of  the  absorbing  gases  present,  Pg,  accounts  for  this  factor,  extension 
to  the  consideration  of  flame  emission  concerns  the  replacement  of  the 
gas  values  for  those  values  found  in  the  flame. 

i 


4  °SB  Tf4  PG  (j  p-- r)/  (3-12) 

Inherent  in  this  extension  is  the  assumption  that  the  flame  consists  of 
only  carbon  dioxide  and  water  vapor.  Such  an  assumption  is  valid  for 
most  combustion  systems  because  COg  and  are  the  only  species  pre¬ 
sent  in  large  quantities  which  have  sufficiently  wide  absorption  bands 
to  attenuate  the  radiation  (McAdams,  1954).  Thus  the  partial  pressure 
is: 


(3-13) 


Hottel  has  also  found  that  the  value  of  the  differential  of  e^.  is  ap¬ 
proximated  by: 


(3-14) 


which  may  readily  be  utilized  in  the  two  dimensional  formulation. 

Using  the  equations  developed  above,  the  one  dimensional  analysis 
may  be  extended  to  consider  two-dimensional  exchange  between  all  elements 
of  the  flame  and  an  area  of  the  wall  (Anon.,  1968).  The  total  exchange 
is  a  function  of  the  volume  element  of  the  flame,  the  emission  per  unit 
volume  (Eqn.  3-12),  the  transmittance  (a  measure  of  radiation  attenuation), 
and  a  view  factor  which  will  be  described  later.  The  problem  considered 
and  the  symbols  used  are  illustrated  in  Fig.  8. 
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The  elemental  volume  of  the  flame  is  considered  to  be  the  elemental 
length  times  the  cross  section  of  the  combustion  chamber  at  this 
point.  Physically  this  seems  to  be  a  poor  approximation,  but  under  the 
assumption  of  one  dimensional  flew  and  a  mean  temperature  across  the 
chamber  the  analysis  is  as  accurate  as  possible.  Therefore,  the  volume 
of  the  flame  element  is: 


(3-15) 


where  denotes  the  cross  section  at  that  axial  position. 

The  view  factor,  F^  ,  may  be  defined  as  the  radiant  energy  arriving 
at  an  element  of  the  wall  divided  by  the  total  energy  leaving  the  flame 
element  in  all  directions.  The  view  factor  is  dependent  on  the  specific 
geometry  of  the  system  and  the  determination  is  possible  for  most  geo¬ 
metries  of  interest  (McAdams,  1954,  Sparrow  and  Cess,  1966).  The  de¬ 
termination  of  specific  view  factors  may  be  found  in  numerous  texts 
(McAdams,  1954,  Sparrow  and  Cess,  1966,  Howell  and  Siegel,  1969)  and 
articles  (Hamilton  and  Morgan,  1952)  and  further  explanation  or  deriva¬ 
tion  would  be  necessary  only  if  a  specific  geometry  were  being  considered. 

The  transmittance  is  a  measure  of  the  energy  emitted  from  the  ele¬ 
mental  flame  volume  which  arrives  at  the  wall  element  unattenuated. 
Attenuation  is  the  result  of  absorption  by  the  intervening  gases.  Hottel 
(McAdams,  1954)  has  provided  typical  empirical  curves  depending  on  the 
partial  pressure  of  the  gases  present  and  the  distance  between  the 
element  of  flame  and  the  element  on  the  wall,  s.  Northern  Research 
(Anon.,  1968)  has  correlated  this  result  as  shown  in  Fig.  9  using  the 
followi ng: 


14.82 

PQ  s  +  14.82 


(3-16) 


An  average  or  mean  value  of  the  transmittance  F  ,  has  also  been  de- 
3  mean 

termined  to  account  for  the  variation  of  s  over  the  radius  of  the  com¬ 
bustor  for  an  annular  combustor: 
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for  which 

rl  =  PG  [r3  ’  r2  "  (1/2°  "  1/4°Hrl  ~  r2} 


(3-17) 


where  n  Indicates  the  number  of  elements  across  the  radius  for  which  s 
was  averaged. 

Introducing  the  above  variation  into  Eqn.  (3-5)  accounts  for  the 
emission  of  the  flame  element  to  the  wall: 


1  +  a  de,. 

:~rJ1)pG(jTp:)0 


(3-18) 


The  total  radiant  energy  reaching  the  wall  element  will  be  the  sunmation 
over  all  flame  positions: 


The  above  formulation  thus  accounts  for  radiation  reaching  the  wall 
element  from  all  of  the  flame  positions. 

In  conclusion,  the  methods  of  radiation  transfer  theory  were  not 
applied  to  the  problem  of  flame  radiation  because  of  the  mathematical 
complexity  involved.  Winter  (1955)  has  provided  a  method  to  approxi¬ 
mate  the  radiation  exchange  between  an  element  of  the  flame  and  an  ele¬ 
ment  of  the  wall  at  the  same  axial  location.  Finally,  the  extension 
to  determine  the  transfer  from  the  entire  flame  to  all  axial  positions 
was  made  using  the  analyses  of  Hottel  (McAdams,  1954)  and  Northern 
Research  and  Engineering  Corporation  (Anon.,  1968).  Only  one  factor 
need  be  determined  to  calculate  the  radiation  exchange  using  Eqn.  (3-5), 
the  flame  emissivity.  Determination  of  the  ernissivity  of  both  luminous 
and  nonluminous  flames  is  the  subject  of  the  next  section. 
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C.  Flame  Emisslvlty 

The  determl natl on  of  flame  emlsslvltles  has  been  an  empirical 
matter  In  the  majority  of  studies.  The  most  widely  used  results  are 
those  of  Hottel  (McAdams,  1954)  for  carbon  dioxide  and  water  vapor, 
the  assunptlon  being  that  the  flame  consists  mainly  of  these  two  emit¬ 
ters.  This  Is  true  only  for  the  case  of  non! uml nous  flames.  Luminosity 
Is  the  result  of  the  presence  of  solid  carbon  particles  emitting  con¬ 
tinuous  radiation.  The  dlstlnltlon  between  the  two  types  of  flames, 
the  approximations  made  In  the  literature  to  calculate  the  enrfsslvltles 
of  both  categories  from  empirical  results,  and  a  brief  review  of  theo¬ 
retical  approaches  to  the  problem  are  discussed  below. 

The  empirical  charts  prepared  by  Hottel  (McAdams,  1954)  are  ad¬ 
vantageous  because  they  account  for  the  emission  of  real  gases  (Einstein, 
1963).  Calculation  of  the  emissivlty  from  Hottel 's  results  depends  on 
the  partial  pressure  of  C02  and  H20  present  In  the  flame,  the  charact¬ 
eristic  radius  of  the  geometry  considered  L,  the  flame  temperature, 
and  the  wall  temperature.  The  necessity  of  providing  both  the  flame 
temperature  and  the  wall  temperature  Is  based  on  Klrchcff's  1  va  (Winter, 
1955);  that  Is,  the  absorption  of  energy  by  the  gas  depends  on  the  tem¬ 
perature  of  the  emitter  (In  this  case  the  wall  at  T  )  and  emits  energy 

n 

at  the  temperature  of  the  gas,  Tq.  Correction  factors  have  been  pro¬ 
vided  to  account  for  changes  In  pressure  and  the  variation  In  mole 
fractions  of  water  and  carbon  dioxide  present.  One  disadvantage  of  the 
results  presented  is  that  they  are  for  one  atmosphere  and  extreme  ex¬ 
trapolation  Is  required  to  obtain  emissivities  and  absorptivities  for 
typical  gas  turbine  pressures  (Anon.,  1968).  A  second  disadvantage  Is 
that  for  typical  fuels  and  local  fuel  to  air  ratios,  luminous  carbon 
emission  must  somehow  be  added  to  the  basic  nonluminous  radiation.  To 
provide  results  for  various  partial  pressures  of  C02  and  H2C  and  empir¬ 
ical  results  over  a  range  of  carbon  particle  concentrations  would  be 
a  difficult  task.  For  this  reason,  most  Investigations  have  chosen  to 
use  continuous  functions  to  predict  flame  properties  (Lefebvre  and  Her¬ 
bert,  1960,  Tlpler,  1955,  Anon.,  1968,  Aref  and  Sakkal ,  1962,  Thri ng 
et  al.,  1961,  Beer  and  Claus,  1962). 

A  contlnous  functional  relationship  has  been  found  to  correspond 
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to  Hottel's  charts  within  one  or  two  per  cent  (Anon.,  1968).  For  non¬ 
luminous  flames  the  parameters  on  which  the  value  of  a f  and  depend 
were  included  in  relations  suggested  by  Lefebvre  and  Herbert  (1960): 

ctf  -  (rl)  *2  (P)  <f>3  (Tw)  (3-20) 

ef  ~  4>1  (rl)  <j>2  (P)  <f>3  (Tf)  (3-21) 

Each  of  these  factors  were  incorporated  in  the  relation  used  by 
Lefebvre  and  Herbert  (1960)  which  was  provided  by  Reeves  (Lefebvre 
and  Herbert,  1960). 

m  0.4  -1.5 

ef  =  1  -exp  [1.60  x  104  P  (r  if)  Tf  ]  (3-22) 

For  the  remainder  of  the  discussion  the  expression  for  absorptivity  will 
be  omitted;  the  only  change  necessary  is  to  replace  the  flame  tempera¬ 
ture  with  the  source  temperature,  T  .  However,  Eqn.  (3-22)  only  pro- 
vides  results  for  the  nonluminous  flame  and  not  the  luminous  results 
which  are  of  interest  in  some  gas  turbine  studies. 

The  distinction  between  the  emissions  on  a  spectral  basis  of  a 
luminous  and  nonluminous  flame  is  illustrated  in  Fig.  10  and  11  from 
the  work  of  Weeks  and  Saunders  (1958).  The  emission  rate  variation  over 
the  spectrum  substantiates  the  assumption  made  about  the  flame  consisting 
primarily  of  H20  and  CO2,  excluding  carbon  monoxide.  The  emission  bands 
of  carbon  dioxide  are  found  at  1 . 9p  and  2.7y.  The  2.7p  band  is  over- 
laped  with  the  H2O  emission  band  which  has  several  peaks  in  the  2.7y 
range.  The  presence  of  water  is  identified  by  the  peak  just  below  2.6y. 

If  carbon  monoxide  were  present  in  sufficient  concentrations  to  attenuate 
radiation,  there  would  be  a  peak  in  the  2.3  to  2.4u  range.  Thus,  it  is 
concluded  that  for  a  typical  fuel  such  as  kerosene  the  omission  of  CO 
from  consideration  is  acceptable. 

Fig.  11  shows  the  spectral  variation  in  emission  rate  due  to  the 
presence  of  carbon  particles  for  a  typical  heavy  composite  fuel  ("Mothball'’). 
Two  important  factors  are  illustrated  in  this  figure:  firstly,  radia¬ 
tion  from  the  solid  carbon  particles  is  spectrally  continuous  rather 
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than  In  discrete  bands  as  found  for  gas  emissions;  secondly,  the  emission 
rate  of  carbon  particles  Is  considerably  greater  than  that  of  the  non- 
lumlnous  HgO  and  C02  radiation.  The  peak  illustrated  in  Fig.  10  for  C02 
and  H20  in  the  2.7u  range  Is  seen  on  the  right  hand  side  of  Fig.  11. 

The  variation  In  magnitude  Is  due  to  the  higher  heat  release  value  of 
mothball  fuel  than  for  kerosene.  In  the  discussion  to  follow.  Fig.  11 
Illustrates  why  some  authors  have  chosen  to  neglect  the  nonluminous  radi¬ 
ation  entirely. 

Lefebvre  and  Herbert  (1960)  were  the  first  to  account  for  the  vari¬ 
ation  In  emlsslvlty  due  to  the  luminous  components  with  the  addition  of 
a  constant  factor  to  the  emlssivity  equation,  (Eqn.  (3-22)). 

ef  =  1  -exp  [1.6  x  104  A  P(rF)°*5  Tf'1*5]  (3-23) 

The  factor  A  is  strictly  an  empirically  determined  factor.  Typical 
values  for  various  fuels  are  shown  In  Table  2  (Lefebvre  and  Herbert, 
1960).  Note  the  difference  In  the  luminous  factor  for  the  two  fuels, 
kerosene  and  mothball,  shown  in  Figs.  10  and  11,  Lefebvre  and  Herbert 
determined  that  the  luminosity  factor  A  could  be  correlated  to  the 
carbon  to  hydrogen  ratio  of  the  fuel  In  use  by  the  following; 

A  »  7.53  (C/H  -  5.5)0*84  (3-24) 

Northern  Research  (Anon,,  1968)  has  attempted  to  correlate  the  luminosity 
factor  to  provide  results  more  accurate  than  Lefebvre  and  Herbert's. 
Firstly,  an  exponential  relation  of  the  form: 

A  =  exp  (C/-gj-‘4)  (3-25) 

was  tried  and  the  results  are  shown  in  Fig.  12  compared  to  the  data  of 
Schirmer  and  Quigg  (Anon.,  1968).  A  recorrelation  of  the  results  has 
been  provided  of  the  form: 

ef  »  1  -exp  [4.7  p1*3  A  (rl)0,5  Tf"1,5] 


(3-26) 
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Table  2. 

Luminosity  Factors  for  Various  Fuels 
(Lefebvre  and  Herbert,  1960) 


Fuel 

A 

Kerosene 

1.7 

Gas  Oil 

6.7 

Mothbal 1 

16.2 

Tarmac 

19.1 

tamer*  iimiiins 
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amiMcarii.  data  (aw»..  i«b) 

Cat*  fro*  Sdrinaar  *as  Qoigg  (A am..  1968}  wuaiig  Tft  -  3600  <ji*g 
Valves  predicted  fruo  Is t&C  1966  correlation  (fgn.  3-?t  ) 

Values  predicted  f rat  ttifaC  1%4  correlafoo  {Cq*.  >») 
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where 


A 


0.74 


The  results  for  the  revised  correlation  which  uses  a  luninosity  factor 
similar  to  the  Lefebvre  correlation  with  a  power  law  dependence  on 
pressure  in  the  exponential  term  are  also  presented  Fig.  12. 

In  an  attempt  to  provide  a  physical  basis  for  the  luminosity  factor, 
Hottel  (McAdams,  1954),  Thring  et  al .  (1962),  and  BeeV  and  Claus  (1962), 
have  suggested  expressions  of  the  form: 


cf  ■  1  -exp  t-C  l  f  (X)]  (3-27) 

where  C  is  the  soot  or  carbon  particle  concentration,  i  the  mean  free  path, 
and  f  (x)  Is  a  function  of  the  wavelength,  particle  size,  and  particle 
distribution.  Thring  et  al.  (1962)  assumed  that  the  luninous  flame  be¬ 
haved  similarly  to  a  turbid  medium.  By  application  of  light  scattering 
theories,  the  authors  claim  to  have  derived  the  following  relation: 


cf  =  1  -exp  [-Kj  Cl]  (1  -  e  )  (3-28) 

where  C  is  the  average  soot  concentration.  This  is  in  contrast  to  the 
formulation  found  from  empirical  results: 


ef  =  1  'exp  £~K1  Cmax  ^  ^  '  S'  (3"29^ 

A 

where  C|J)ax  is  the  maximum  soot  concentration  and  accounts  for  the 
nonluminous  radiation  factor.  Thring  et  al .  assumed  the  particles  were 
in  one  of  two  classes:  spherical  units  of  from  100  to  800  A  with  a 
nominal  size  of  400  A  or  spheres  from  1  to  15y, 

Stull  and  Plass  (1960)  investigated  the  interaction  of  electromagnetic 

O 

waves  with  particles  ranging  from  50  to  1000  A,  developing  the  optical 
theories  by  considering  the  particles  to  be  equally  dispersed.  The  re¬ 
stricting  assumptions  made  by  Sato  et  al .  (1962)  are  typical  of  th  s 
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type  of  investigation.  The  particles  were  assumed  to  be  identical  i rs 
size  and  material  and  to  be  equally  dispersed  throughout  the  flame.  The 
temperature  of  the  cloud  was  also  assumed  to  be  uniform.  The  physical 
Insight  gained  by  investigating  such  a  system  rather  than  assuming  a 
constant  factor  of  luminosity  is  lost  in  the  idealities  of  the  model 
proposed.  This  criticism  is  specifically  pointed  to  the  assumption  of 
a  uniform  dispersion  and  temperature  In  the  flame. 

In  conclusion,  with  the  present  knowledge  concerning  carbon  forma¬ 
tion  in  flames,  the  assumptions  necessary  to  arrive  at  an  entirely  an¬ 
alytical  solution  have  little  practical  applicability.  The  accuracy 
and  convenience  of  the  continuous  function  appear  to  have  overshadowed 
the  empirical  charts  of  Hottel,  especially  with  the  addition  of  luminous 
radiation.  However,  the  expressions  used  must  be  verified  as  to  the 
accuracy  attainable  when  extrapolated  to  higher  pressures. 

D.  Radiation  Exchange  Between  the  Chamber  Halls 

The  consideration  of  radiation  exchange  between  different  areas  of 
the  wall  was  originated  in  the  work  of  Tipler  (1955).  The  specific 
case  considered  was  the  exchange  between  the  wall  in  the  primary  zone 
and  that  of  the  secondary  zone.  A  similar  analysis  was  derived  to  de¬ 
termine  the  general  exchange  between  wall  areas  for  an  annular  combustor 
(Anon.,  1968). 

The  analysis  as  proposed  by  Tipler  suggests  that  the  exchange 
between  the  primary  wall  of  a  tubular  combustor  and  the  secondary  wall 
should  be  given  by  an  equation  similar  to  Eqn.  (3-5): 

4  R3  -  “SB  FdA,-dA2  tS,  Tw,4  -  \ 


where  a,.  and  e,  will  be  approximate  the  same  and  the  view  factor, 
^dAi-dA2  depend  on  the  geometry  and  type  of  the  conbustor.  Appro- 

-pri ate  View  factors  are  found  in  numerous  heat  transfer  texts  (Sparrow 
and  Cess,  1966,  McAdams,  1954).  Tipler  included  an  additional  factor, 

(1  +  Cf),  which  accounts  for  the  absorption  by  the  flame.  The  situation 
considered  is  illustrated  in  Fig.  13  and  follows  the  relation: 


<'  -  “f>  FdA,-dA2  “SB  Tw/  -  «w  Tw241 


(3-31) 


M 


FIGURE  13  RADIATION  INTERCHANGE  BETWEEN 
FLAME-TUBE  WALLS 
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Since  a,  =  -  for  most  conditions: 

w  w 


^3  s  ^  ‘  cf*  cw  aSB  FdA,-dA0  £Tw.4  " 

;  u  1  c 


\o- jc  ; 


Northern  Research  (Anon.,  1968)  has  provided  a  similar  analysis  which, 
rather  than  accounting  for  a  flame,  accounts  for  the  absorption  of  the 
intervening  gas  between  the  two  areas.  Under  the  assumption  that  the 
flame  consists  of  C02  and  HgQ  the  two  analyses  are  conceptually  the 
same.  The  transmittance  of  radiation,  r,  is  determined  with  Eqn.  (3-16). 


AR~  =  Ocr  e.  F,,  r  (T  4  -  T  4) 

3  SB  w  dA^-dAg  Wg 

Summation  over  all  wall  positions  yields: 


0SB  ew 


Z  F 

All  Wall  “"run2 
Positions 


dA.-dA0  r  (Tw,  “  Tw„  ^ 
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(3-33) 


(3-34) 


Similar  expressions  represent  the  exchange  between  each  of  the  wall 
elements  and  the  remainder  of  the  system. 

E.  Radiation  Exchange  Between  the  Chamber  Wall 
and  the  Casing 

With  the  basic  principles  of  radiation  exchange  having  been  pre¬ 
sented,  the  presentation  of  the  methods  used  to  determine  R^  is  re¬ 
latively  brief.  The  only  variation  found  in  the  formulation  is  the 
replacement  of  the  wall  temperature  in  Eqn.  (3-30)  with  the  casing 
temperature  tca  (Anon.,  1968,  Lefebvre  and  Herbert,  1960). 

R2  "  FW*C  °S8  <Tw4  “  TCA4)  (3'35) 


Two  assumptions  are  not  Immediatly  evident  in  Eqn.  (3-35).  First,  the 
view  factor  must  account  for  direct  exchange  and  exchange  through  re¬ 
flection  and  rereflection.  A  typical  expression  has  been  given  by 
Northern  Research  (Anon.,  1968)  as: 


wc 
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(3-36) 
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where : 

F7ZT  2  hi ack- surface  overall  view  factor 

MU 

A  ,  h  area  of  the  wall 
w 

Ara  =  area  of  the  casing 

eq  s  emissivity  of  the  casing  (eCA  =  aCA) 

Second,  the  gases  present  In  the  annulus  do  not  attenuate  radiation. 
This  corresponds  to  the  fact  that  02  and  N2  do  not  absorb  or  emit 
radiation  in  the  infrared  (McAdams,  1954).  The  summation  necessary  for 
the  chamber  wall  interchange  is  not  necessary  because  the  casing  is 
asswned  to  have  a  constant  temperature  (Anon.,  1968).  The  variation 
due  to  combustor  geometry  affects  only  the  exchange  factor. 
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CHAPTER  IV 

CONVECTIVE  HEAT  TRANSFER  IN  COMBUSTION  CHAMBERS 
A.  Introduction 

The  evaluation  of  the  heat  transfer  terms  C^  and  C2  has  received 
considerably  less  attention  than  the  problem  of  radiation  in  the  com¬ 
bustion  chamber.  Methods  for  calculating  the  internal  convection  near 
the  front  of  the  chamber  have  been  presented  in  the  literature,  but  the 
impression  left  is  that  the  particular  Investigators  realized  convective 
heating  occurred  and  thus  a  term  accounting  for  it  was  necessary.  The 
problem  centers  around  the  diversity  of  recirculation  flows  found  in 
the  region  and  the  prediction  of  their  magnitude  and  direction.  In 
sharp  contrast  is  the  vast  number  of  studies  concerning  film  cooling, 
a  variation  of  internal  convection.  Unfortunatly  most  studies  have 
been  performed  with  either  low  mainstream  temperatures,  on  the  order  of 
100  to  200°C,  or  with  heated  slot  flow  and  cold  mainstreams.  The  use 
of  these  correlations  to  predict  extrapolated  conditions  could  be  ex¬ 
tremely  misleading.  Spalding  (1965)  has  provided  a  semi -empirical  theory 
which  indicates  the  correct  trends  and  provides  a  basis  of  comparison 
with  experimental  results.  The  final  source  of  internal  convection, 
the  penetration  jet,  has  not  been  considered  at  all  in  any  of  the  heat 
transfer  analyses  found,  for  reasons  which  will  be  discussed  later. 

The  first  portion  of  this  chapter  concerns  the  review  of  the  literature 
written  before  1966  regarding  these  Internal  modes  of  convection. 

External  to  the  flame  tube,  the  convective  cooling  C2  is  a  much 
easier  mode  to  analyze*  The  reduction  in  flow  velocities  are  relatively 
easy  to  determine  and  the  cooling  is  not  dissimilar  to  flow  over  a  flat 
plate.  The  difference  in  the  difficulty  In  determining  internal  and 
external  convective  contributions  is  the  difference  found  in  modelling 
the  flow  conditions  in  the  two  regions.  The  final  topic  to  be  covered 
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Is  thus  the  analysis  of  the  convection  due  to  external  flow,  C0. 

<u 

E.  Internal  Convection,  C-j 

B.  1)  Convective  Heating  in  the  Absence  of  Secondary  Flow 

In  the  front  portion  of  the  primary  zone  is  a  region  in  which  little 
secondary  air  has  been  added  through  penetration  jets  or  film  cooling 
slots.  The  gases  present  have  either  entered  through  holes  in  the  done 
or  have  recirculated  in  a  manner  similar  to  that  shown  in  Fig.  2.  The 
difficulties  encountered  in  modelling  this  flow  have  been  discussed  pre¬ 
viously  (Hammond  and  fie  11  or ,  1971).  Of  the  analyses  concerning  heat 
transfer  in  combustion  chambers  (Lefebvre  and  Herbert,  1960,  Tipler,  1955, 
Winter,  1955),  only  Lefebvre  and  Herbert  have  even  considered  the  effects 
of  flow  in  this  region.  Due  to  the  stratified  flow  in  the  primary  zone 
and  the  fact  that  gases  are  undergoing  rapid  physical  changes  in  this 
region,  the  authors  admitted  that  an  account  of  the  convection  would  at 
best  be  highly  approximate.  The  analogy  proposed  is  to  model  the  flow 
as  straight  pipe  flow.  The  expression  for  such  a  flow  is: 

no  Qn  d*  -0.2 

St  =  0.0283  Re  *  0.0283  (jf- -f)  (4-1) 

or 

Q  0.8  (T  -  T  ) 

C,  =  0.020  a  o.8  (t&)  -4r-a-y-  (4-2) 

M  p  Up 

where  Qp  is  the  mass  flow  rate,  Ap  the  cross  sectional  area,  and  d* 
the  hydraulic  diameter,  all  of  the  primary  zone.  Immediately  upon 
statement  of  the  above  formulation  Lefebvre  and  Herbert  express  some 
doubt  as  to  the  validity  of  the  equation.  Firstly,  the  flow  velocity 
near  the  wall  is  found  to  be  higher  than  the  mean  velocity  in  the  pri¬ 
mary  zone  because  of  the  necessity  of  such  a  flow  to  stabilize  the 
flame.  Secondly,  the  use  of  swirl  further  increases  the  approxinate 
nature  of  the  solution.  Finally,  the  temperature  of  the  gases  near 
the  wall  is  lower  which  will  lead  to  higher  gas  densities  near  the  wall 
and  further  exaggerate  the  velocity  effects  mentioned  before. 

An  additional  problem  is  found  in  determining  the  temperature  at 
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which  the  properties  of  the  fluid  should  be  analyzed.  For  radiation, 
the  bulk  gas  temperature  may  be  satisfactory,  but  due  to  the  stratifica¬ 
tion  of  temperature  in  the  flow,  the  value  of  Tr  for  convection  will 

u 

be  somewhere  between  and  T  . 

Calculation  of  the  convective  heating  in  the  above  manner  after 
secondary  air  has  been  added  will  becone  a  progressively  better  approxi¬ 
mation  at  positions  down  the  chamber  since  the  recirculation  flow  will 
disappear  and  the  flow  become  more  analogous  to  pipe  flow. 

In  conclusion,  the  determination  of  in  the  primary  zone  has  re¬ 
ceived  very  little  attention  in  the  literature.  The  flow  has  been 
modelled  as  being  similar  to  pipe  flow.  The  error  involved  in  such  an 
analogy  will  depend  on  how  closely  pipe  flow  resembles  the  true  situation 
in  the  combustor  geometry  considered.  Fortunately,  a  large  portion  of 
the  secondary  zone  region  is  cooled  by  film  and  splash  cooling  devices 
for  which  the  cooling  effects  are  better  understood. 

B.  2)  Film  Cooling  of  Chamber  Walls 

To  protect  the  chamber  wall  from  convective  heating,  design  practices 
have  dictated  the  use  of  film  cooling  slots.  The  film  cooling  slot 
(Fig.  14)  provides  a  blanket  of  cool  annulus  air  to  protect  the  liner 
wall  from  the  hot  chamber  gas.  The  cooling  flow  may  consist  of  three 
distinct  zones:  the  potential  core,  the  wall  jet  region,  and  the  fully 
developed  turbulent  boundary  layer  region.  The  character  and  the  ex¬ 
istence  of  these  regions  depends  on  the  ratio  of  coolant  or  slot  gas 
velocity,  u$,  to  the  magnitude  of  the  free  stream  velocity,  um$.  For 
slot  velocities  much  greater  than  the  mainstream  velocity,  us  >  >  ums, 
the  wall  jet  velocity  may  be  calculated  utilizing  Glauert's  (1956)  wall 
jet  theory.  For  the  case  where  u  >>  u,.  the  wall  jet  region  no  longer 

no  a 

exists  (Stollery  and  El-Ehwapy,  1965).  The  flow  transition  is  thus  from 
the  potential  core  flow  directly  to  the  fully  developed  turbulent 
boundary  layer  flow.  Because  both  of  the  situations  are  found  in  com¬ 
bustion  chambers,  the  theories  and  experimental  evidence  pertaining  to 
each  will  be  reviewed. 

Of  the  numerous  works  related  to  film  cooling,  the  majority  concern 


FIGURE  14  FILM  COOLING  SLOT  AND  FLOW  REGIONS 
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the  turbulent  boundary  layer  region;  considerably  less  work  has  been 
done  on  the  wall  jet  flow  and  Aimoet  the  potential  cote.  In 

the  case  of  the  v/all  jet,  the  theory  of  Glauert  (1956)  will  be  reviewed 
and  the  success  of  the  theory  compared  to  the  experimental  findings  of 
Seban  and  Back  (1960).  Spalding  (1965)  has  proposed  a  semi -empirical 
theory  which  encompasses  both  the  wall  jet  and  turbulent  boundary  layer 
region.  The  basis  for  the  theory  will  be  reviewed  and  the  wall  jet 
solution  contrasted  to  that  of  Glauert.  The  solution  for  the  turbulent 
boundary  layer  will  be  stated  as  a  guide  to  presenting  the  numerous 
experimental  investigations  of  the  turbulent  region.  Hopefully,  the 
order  chosen  for  the  review  will  provide  substantially  greater  insight 
Into  the  numerous  correlations  presented  In  the  literature. 

Throughout  the  literature,  the  investigations  (both  theoretical 
and  experimental)  employ  the  "adiabatic  wall"  for  their  studies.  This 
Is  merely  a  reference  state  which  avoids  the  introduction  of  extraneous 
phenomena  into  the  results.  Without  such  a  reference  state,  the  results 
of  the  various  authors  could  not  accurately  be  compared.  Additionally, 
the  scatter  found  in  the  results  of  the  various  investigators  may  be 
partially  attributed  to  insufficiently  adiabatic  walls. 

B.  2-1)  The  Potential  Core  Region 

The  potential  core  exists  from  the  slot  exit  until  the  turbulence 
Induced  by  the  mainstream  mixes  the  slot  and  mainstream  flows  (see 
Fig,  14).  The  gases  below  the  mixing  line  retain  the  properties  of  the 
coolant  flow.  For  this  reason,  in  the  absence  of  radiative  heating, 
the  adiabatic  wall  temperature  will  be  that  of  the  coolant  flow  (the 
effect  of  radiative  transport  is  included  through  the  use  of  the  theory 
of  superpositing  as  described  in  Chapter  III).  Thus,  the  convective 
heat  transfer  is  determined  using  Newton's  convective  heat  transfer 
equation  (McAdams,  1954); 


■  hcA<Tw-V 


(4-3) 


where  the  value  of  hc,  the  convective  heat  transfer  coefficient,  is  that 
of  laminar  flow  over  a  flat  plate,  T  is  the  adiabatic  wall  temperature, 
and  T$  the  temperature  of  the  slot  fluid.  Assuming  the  region  behind 
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the  slot  is  similar  to  a  flat  plate  is  acceptable  because  the  normal 
slot  height  is  surncientiy  small,  compared  to  the  radius  of  thr  chamber, 
that  effects  due  to  the  curvature  will  be  negligible.  An  additional 
assumption  is  that  the  slot  is  continuous  around  the  chamber  to  avoid 
the  Inclusion  of  end  effects.  In  conclusion,  the  potential  core  in  the 
absense  of  radiation  is  unaffected  by  the  presence  of  the  mainstream 
thus  protecting  the  wall  from  convective  heating. 

B.  2-2)  The  Wall  Jet  Region 

The  wall  jet  flow  region  velocity  profile  (Fig.  14)  Indicates  why 
the  wall  jet  must  be  considered  as  a  specific  type  of  flow  rather  than 
as  a  fully  developed  turbulent  boundary  layer.  The  flow  near  the  wall 
is  characteristic  of  the  viscous  effects  of  the  boundary  layer.  The 
unique  property  of  the  flow  Is  the  peaked  velocity  profile  near  ■‘‘he  wall 
which  slowly  dissipates  to  the  freestream  value  as  the  end  of  tin.  velocity 
boundary  layer  is  approached.  Of  course,  such  a  region  does  not  exist 
when  the  coolant  flow  velocity  is  less  than  or  of  the  order  of  the 
magnitude  of  the  mainstream  velocity. 

The  pioneering  study  of  the  wall  jet  was  conducted  by  Glauert  (1956), 
who  solved  the  boundary  layer  problem  for  the  velocity  profile  assuming 
there  is  a  similarity  solution  such  that  the  velocity  u  -  xa  and  the 
velocity  boundary  layer  thickness  5  is  proportional  to  xb,  6  -  xb. 

The  value  of  a  and  b  must  be  found.  The  equations  of  interest  for  the 
system  are  the  two  dimensional  continuity  equation: 


3u  ,  3v 
■Sx  dy 


0 


and  the  momentum  equation: 


u  ax  +  v  ay 


au  „ 


au> 
t  W 


(4-4) 


(4-5) 


where  the  bars  above  the  velocity  components  represent  the  mean  velocity 
of  the  flow.  It  should  be  noted  that  the  eddy  viscosity  has  been  used 
in  place  of  the  two  terms  for  viscosity  in  Eqn.  (2-2).  The  boundary 
conditions  available  for  Eqn.  (4-4)  and  (4-5)  are: 


p  * 

9J 


u  =  v  =  0  at  y  «  0  and  u  u^as  y  +  ®  (4-6) 

For  the  use  of  Glauert's  formulation  it  will  be  required  that  the  main¬ 
stream  velocILy  approacn  2ero. 

To  simplify  the  task,  Glauert  used  Prandtl’s  hypothesis  to  describe 
the  behavior  of  e^.,  that  the  eddy  viscosity  Is  constant  across  the 
boundary  layer  and  proportional  to  the  product  of  the  maximum  mean  velo¬ 
city,  and  6,  the  boundary  layer  thickness. 

Assuming  there  Is  a  similarity  solution  for  the  velocity  profile 
r.r.  u  -  xa  and  <S  -  x*5,  Glauert  determined  that  e  -  xa+\  where  b  * 

1  iuvi  a  -  -4/3.  By  defining  stream  functions  and  integrating  to  solve 
cqr .  '4-5*  '’lauert  found  his  results  gave  the  characteristic  results 
of  tiie  outer  portion  of  the  laminar  wall  jet  shown  In  Fig.  15.  To 
account  for  the  remainder  of  the  flow,  Glauert  Incorporated  the  frictional 
effects  found  in  turbulent  pipe  flow  studies  by  Blasius: 

t  *  0.0225  p  u2  (±jp-)1/4  (4-7) 

that  is,  the  frictional  effect  of  the  wall  has  been  accounted  for  as  a 
function  of  the  parameters  of  the  system.  The  result  of  matching  the 
theoretical  profile  provided  by  Prandtl's  hypothesis  for  the  flow  beyond 
the  maximum  velocity  with  the  empirical  results  of  Blasius  is  shown  in 
Fig.  16,  where  and  f'  are  the  nondimensional  distance  and  velocity 
functions. 

Glauert  also  defined  a  function,  ayj,  which  is  a  measure  of  the 
extent  of  the  fractional  effects  on  the  flow  for  the  velocity  profile 
beyond  the  maximum  velocity.  This  function  will  be  of  use  later  in 
studying  the  empirical  results  of  Seban  and  Back  (1960). 

The  above  is  a  cursory  review  of  the  work  of  Glauert  and  in  no 
way  implies  the  solution  for  the  boundary  layer  is  a  simple  one.  How¬ 
ever,  the  analytical  aspects  have  been  avoided  because  the  experimental 
results  will  provide  needed  information.  The  temperature  profiles  can¬ 
not  be  obtained  from  Glauert’s  numerical  solution  by  a  technique  such 
as  the  application  of  Reynold’s  analogy  (Schlichting,  1960)  because  the 
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FIGURE  15  VARIATIONS  OF  VELOCITY  (f * )  WITH  DISTANCE  FROM  THE 
WALL  (n)  FOR  A  LAMINAR  WALL  JET  (GLAUERT,  1956) 
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flow  and  temperature  fields  are  not  analogous.  Therefore,  Seban  and 
Back's  (1960)  results  have  to  be  used  rather  than  extant  ««  Gl’ucrt's 
theory  to  temperature  profiles. 

These  investigators  empirically  determined  the  velocity  and  temper¬ 
ature  profiles  in  the  wall  jet.  The  analysis  used  strictly  follows  the 
work  of  Glauert  reviewed  above.  The  system  Seban  and  Back  studied  is 
depicted  in  Fig.  17.  The  authors  indicated  that  the  difficulty  in 
applying  Glauert's  results  stems  from  the  velocity  expression: 


where  U  is  a  reference  velocity  which  is  not  easily  defined,  a 
function  which  depends  on  the  ratio  of  the  nondimensional  distance  from 
the  wall  to  the  maximum  velocity,  to  the  nondirajnsional  distance  from 
the  maximum  velocity  to  the  end  of  the  boundary  layer.  The  proportionality 
is  defined  by  which  in  turn  yields  the  power  of  u  -  xd  and  6  ,  xb  by: 


*  +  aWJ 

k  4  +  Hj 


(4-9) 

(4-10) 


or  approximately: 


a  =  UWJ 
V™  "a 


WJ 


(4-11) 

(4-12) 


The  function  df/dn6  is  related  to  the  similarity  variable  n6  which  is 
given  by: 


_  _  4  -  4b  yU  ,Ux> 

hx  -  f— -  fr  (rr) 


-b 


'WJ 


u  'y 


(4-13) 


In  actuality,  the  coordinate  f  =  df/dn6  of  Fig.  16  is  nondimensionalized 
as  follows: 

/(r~) 

V* 
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2ne 


f 


(4-14) 


with  6/2  being  the  distance  from  the  wall  to  where  IT  =  1/2  uj^.  A 
value  of  6/2  is  used  for  the  related  variable  because  this  is  the  easiest 
experimental  point  to  determine  with  sufficient  accuracy. 

Fig.  16  compares  extremely  well  with  the  results  of  Seban  and  Back. 
Generally,  It  was  found  that  as  the  slot  height  was  reduced  the  agreement 
with  the  theoretical  results  Improved.  This  does  not  imply  that  the  re¬ 
sults  were  drastically  far  from  the  theoretical  predictions  at  large 
slot  heights,  becuase  the  accuracy  was  still  fairly  good.  Of  interest 
to  the  application  of  modelling  film  cooling  flow  was  the  fact  that  the 
Influence  of  the  free  stream  flow  was  noted  even  for  mainstream  velocity 
ratios  un/Umax  of  0.40. 

Seban  and  Back  have  correlated  their  data  In  terms  of  6,  which  is 
the  distance  downstream  of  the  slot,  x,  plus  the  starting  or  effective 
starting  length  of  the  boundary  layer,  xQ: 


E  »  x  +  xQ 


(4-15) 


From  Qauert's  results  the  following  analytical  function  was  compared 
with  the  empirical  results: 


xo  +x 


*  6  IT  1/4 

/ o  \  /  max> 

(F™  }  x  * 


4  + 


a, 


[rairc  (1  "  rror“>(aw  +  °-07)  ^6 

f  5/4 

(JS*)  3 

max 


-5,4 


(4-16) 


for  which  fmax,  f^x  and  are  functions  of  a^.  Experimental  results 
Indicate  that  the  power  law  behavior  may  not  exist  far  downstream  of 
the  slot.  This  region  is  of  no  importance  for  practical  film  cooling 
results  and  will  not  be  discussed. 

For  the  experimental  values,  it  was  found  that: 
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_  0.60 

0.212  (rl) 

•rr  • 


(  4-1  7  \ 

,  .  . .  / 


which  agrees  with  experimental  results  within  5  per  cent.  Similarly, 


u  -0,45 

*  3.6  (£r)  (4-18) 

s  s 

for  all  slots  within  a  5  per  cent  error. 

Seban  and  Back  extended  Glauert’s  work  by  analytically  determining 
the  expected  shape  of  the  temperature  profiles.  The  result  Is  compared 
with  the  experimental  profiles  in  Fig.  18  for  a  0.25  inch  slot.  The 
profiles  near  the  slot  approach  a  zero  slope  near  the  wall  as  would  be 
expected . 

0.  B.  Spalding  (1965)  developed  a  theory  which  Includes  all  possible 
Interactiois  of  film  cooling  slots  with  a  flat  plate.  In  actuality, 
the  theory  was  developed  to  include  mass  transfer  through  a  flat  plate 
and  combustion  in  the  boundary  layer,  as  well  as  film  cooling  flow. 

Because  of  Its  substantial  length,  the  theory  will  only  be  briefly  out¬ 
lined:  the  results  for  the  wall  jet  will  be  presented  as  the  final  part 
of  this  section  and  the  results  for  the  turbulent  boundary  layer  will 
lead  Into  tne  subsequent  discussion  of  film  cooling. 

The  "Unified  Theory"  of  Spalding  was  an  attempt  to  organize  the  ex¬ 
tensive  experimental  work  done  on  film  cooling  (Chin  et  al.,  1958,  Chin 
et  al.,  1961,  Hartnett  et  al.,  1961,  Hatch  and  Papell,  1961,  Papell,  1960, 
Papell  and  Trout,  1959,  Seban,  1960,  Seban  and  Back,  1962,  Back  and  Seban, 
1962,  Seban  et  al.,  1957,  Seban  and  Back,  1960,  Stollery  and  El-Ehwany, 
1965).  Each  of  these  investigators  have  varied  the  mainstream  velocity 
or  temperature,  the  coolant  flow  and  temperature,  correlated  the  results, 
and  presented  this  correlation  usually  with  use  of  very  little  physical 
insight.  A  second  difficulty  with  the  experimental  results  is  that  the 
conditions  under  which  the  film  cooling  was  studied  were  usually  at  low 
mainstream  temperatures  (less  than  700°R)  or  with  a  cold  mainstream  and 
heated  slot  flow.  Therefore,  the  correlations  determined  cannot  be  re¬ 
liably  extended  to  the  range  of  interest  in  combustors. 

Fig.  19  is  a  graphical  representation  of  the  possible  situations 
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FIGURE  19  PROCESSES  CONSIDERED  BY  SPALDING  (1965) 
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which  Spalding  has  tried  to  account  for  In  his  theory.  The  purpose  of 
the  unified  theory  was  to  develop  a  set  of  equations  such  that  each  term 
or  set  of  terms  account  for  one  of  the  flow  characteristics  shown.  The 
engineer  could  then  eliminate  those  terms  which  are  not  necessary  for  the 
particular  case  he  is  considering.  In  this  Instance,  all  but  those  ter— s 
pertaining  to  film  coolant  flow  may  be  eliminated  from  consideration. 

The  theory  Is  based  on  two  postulates.  First,  the  profiles  of  velo¬ 
city,  temperature,  and  concentration  are  described  by  expressions  with 
two  main  components,  a  boundary- layer  component  and  a  wall  jet  component. 
The  boundary- layer  component  accounts  for  the  momentum,  heat,  and  mass 
transfer  to  the  wall.  The  jet  component  models  the  effects  of  interaction 
between  the  coolant  stream  and  the  mainstream.  The  second  postulate 
assumes  that  fluid  Is  entrained  Into  the  wall  layer  in  the  same  manner 
as  In  a  turbulent  jet. 

There  are  two  restlrctlons  on  the  following  development:  firstly, 
the  properties  of  fluid  (density,  viscosity,  specific  heat,  and  thermal 
conductivity)  are  constant  across  the  boundary  layer;  secondly,  the  wall 
Is  hydrodynamlcally  smooth. 

The  mathematical  quantities  which  are  of  concern  are  defined  by  three 
Integral  quantities  (I)  and  five  Reynolds  numbers.  The  three  intergral 
properties  depend  on  Z  (=  u/u^)  or  4>  (a  conserved  property)  as  follows: 

I,  '  W»ns>  2  d  *  <4-'9> 

h  5  £  z2d«  (4-20) 

3  t'  2  d  5  (4-2!> 

where  £  =  y/yms,  and  tlie  subscript  ms  stands  for  the  mainstream 
quantity.  These  expressions  may  be  related  to  the  displacement  thick¬ 
ness  6^,  the  momentum  thickness  <52»  and  shape  factor  H: 

<$./y  «  1  -  I, 

rJms  1 


(4-22) 
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W  =  M'23> 

H  *  &}/&2  -  (1  -  I1)/(I1  -  I2)  (4-24) 

In  order  to  provide  conservation  laws  In  the  form  presented  by  Spalding, 
it  Is  necessary  to  define  five  different  Reynolds  nunbers: 


fte  = 

ms 

Dms  ums 

(4-25) 

Re2  = 

p«s  ums  4A«  '  <‘l  ' 

‘  ‘2> 

(4-26) 

Rem  = 
m 

^"S  (cu^S>  ‘  *1 

Rems 

(4-27) 

Remax 

'  pms  “irax  ^max^s 

(4-28) 

Rex  5 

Jf  Kt  VV  dx 

(4-25) 

Using  the  above,  the  conservation  laws  nay  now  be  developed. 

Three  conservation  laws,  conservation  of  mass, momentum,  and  a  con¬ 
served  property  <J>,  may  now  be  written  as  first  order  ordinary  differential 
equations;  for  conservation  of  mass 
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The  conservation  of  raoiiantiitn  equation  may  be  reprpspntPfi  by: 


d  Re-  d  (in  u....  ) 

+  {l  +  H}  Re2_d  -  , 

X  X 


*e244 


d  (*n  w) 

~^r 


«  s  +  m  (4-33) 


where 


^pms  Has  ) 


vf 

T 


(4-34) 


which  Is  a  measure  of  the  shear  force  acting  on  the  fluid. 

The  final  conservation  equation  Is  the  conservation  of  a  property 
such  as  enthalpy  or  species  concentration: 


dlkj  llf  Rem)  +  Rem  "  01  (*T  “  4,ms)  (4"35) 

where  the  subscript  T  stands  for  the  "transfered-substance  state". 

Before  the  three  conservation  equations  may  be  solved*  several 
more  relations  are  necessary.  The  first  of  these  is  the  assumed  velo¬ 
city  distribution  under  uniform-density  conditions: 


Z  *  S1/2  u*  +  (1  +  ZE)  (  -  cos  ir  £)/2 


(4-36) 


The  first  term  of  the  above  relates  to  the  "law  of  the  wall"  and  the 
second  to  the  "law  of  the  wake".  The  terminology  here  relates  to  a 
paper  by  Coles  (Schllchtlng,  1960)  In  which  the  velocity  profile  was 
assumed  to  be  represented  by  a  linear  combination  of  two  universal 
functions,  "the  law  of  the  wall"  being  the  first,  the  "law  of  the  wake" 
the  second.  In  the  above  the  shear  stress  term  function  u+  Is  related 
to  the  nondlmenslonal  distance  from  the  wall  y+,  where: 

y*  s  y{rp)^2/p  (4-37) 


The  quantity  (1  -  Z^)  measures  the  magnitude  of  the  free-mixing -layer 
component  of  the  velocity  profile.  The  final  quantity  (1  -  cos  n  fJ/2 
Is  a  correlation  of  empirical  results. 
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Through  a  series  of  approximations,  Spalding  replaced  the  unknown 
function  u+  and  obtained: 

Z  *  2.5  (s  +  m  ZE)1/2  in  C  +  ZE  +  (1  -  Z£)* 

•(  1  -  cox  If  0/2  (4-38) 


The  Inclusion  of  the  mass  flux  term  through  the  wall  represents  the  shear 
stress  component  due  to  such  a  mass  flux. 

The  second  necessary  relation  Is  the  mass  entrainment  relation, 
which  is  defined  as: 
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ms 


3  +  Z  y 
o  ms 

x 


(4-39) 


where  ZQ  is  related  to  the  velocity  of  the  mixing  layer  and  that  of  the 
free stream. 

Now  that  general  equations  have  been  provided  for  the  physical 
situation  believed  present,  rather  than  continuing  with  the  lengthy 
derivation  of  the  theory.  It  will  be  beniflclal  to  apply  Eqn.  (4-30) 
through  (4-35)  to  the  wall  jet.  The  mass  flux  through  the  wall  is  set 
equal  to  zero  because  the  situation  conslderes  flow  only  from  a  slot. 
The  term  d  (2,n  w)/d  Rex  becomes  zero  if  a  constant  width  of  the  stream 
Is  assumed.  This  reduces  equation  (4-30)  to: 


d  Rem 

OeJ’  =  "  mms  (4_40) 

Similarly,  since  ums  is  assumed  to  be  nearly  zero  for  the  wall  jet,  the 
d  (£n  u^J/d  Rex  term  may  be  eliminated  from  Eqn.  (4-33)  and  (4-35): 
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Using  the  above  relations.  Spalding  showed  that* 
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ID 


P  «c 


■R]7iT 


-  1.837  (C2  At) 
s 


1/2 


(4-43) 


where  the  slot  height  Is  given  by  and  us  represents  the  coolant  velo¬ 
city.  In  the  case  of  the  wall  jet: 


(4-44) 


such  that 


*  v  V 2 

Tp-tp—  «.  1 .837  (Cp  tr)  (4-45) 

s  S  S 

which  using  an  expression  given  by  Spalding  for  6  yields: 


v  -  V2 

*  0.695  (C2  |j£)  (4-46) 

S  5 

Comparison  of  Eqn.  (4-45)  and  (4-46)  to  those  of  Seban  and  Back  (see 
Eqn.  (4-16)  and  (4-17))  indicates  the  theory  is  deflnltly  of  the  correct 
form.  Experimental  results  presented  by  Spalding  show  that  C2  has  been 
determined  by  Stratford  et  al.  so  that  Eqn.  (4-46)  becomes: 


(4-47) 


The  excellent  agreement  exhibited  here  encourages  the  use  of  Spalding's 
enthalpy  flux  equation: 
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(4-48) 


where  fc  =  9.94,  an  empirically  determined  value,  and  n  Is  a  quantity 
which  accounts  for  the  fact  that  heat  and  mass  are  transfered  faster  In 
a  free  jet  than  is  momentum.  Thus,  if  n  =  1,  the  exchange  is  equal. 
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that  is,  heat  is  transferred  at  the  sane  rate  as  momentum.  If  n  =  0, 
heat  and  mass  transfer  is  assumed  to  be  infinitely  fast  compared  to 
momentum  transfer.  Inserting  the  value  of  i.  =  9.94  and  Eqn.  (4-43) 
into  (4-48): 
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(4-49) 


where  h  is  the  enthalpy  of  the  gas  next  to  the  wall  (the  gas  temperature 
can  be  assumed  to  be  very  close  to  the  wall  temperature).  The  value  on 
the  left  hand  side  of  Eqn.  (4-49)  is  a  measure  of  the  efficiency  of  the 
film  coolant  in  maintaining  a  low  wall  temperature.  In  a  latter  part 
of  this  section  it  will  be  more  often  termed  the  efficiency  of  film 
cooling.  Often  the  expression  is  written  as  follows: 
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assuming  the  values  of  Cp  are  relatively  close. 


(4-50) 
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Since  the  above  is  the  result  of  a  semi-empirical  theory,  one  might 
question  the  improvement  over  the  findings  of  experiments  and  correlations 
such  as  those  presented  by  Seban  and  Back  (1960).  The  improvement  lies 
In  the  relatively  few  unknowns  In  Eqn.  (4-51 ).  Rather  than  providing 
a  different  correlation  for  each  mainstream  and  coolant  velocity  ratio, 
only  a  value  of  n  need  be  determined  to  calculate  the  effectiveness  of 
the  film  cooling  wall  jet. 

In  conclusion,  the  work  of  Gleuert  provides  an  understanding  of 
the  flow  issuing  from  a  wall  jet.  However,  since  Reynold's  analogy  may 
not  be  used  to  extend  Glauert's  results,  enpirical  Information  is  neces¬ 
sary  to  ascertain  the  temperature  distribution.  The  theory  of  Spalding 
was  then  developed  in  general  and  applied  to  the  specific  case  of  the 


wall  jet.  The  previously  reviewed  findings  of  Seban  and  Back  (1960) 
substantiate  the  theoretical  predictions  of  Spalding. 

B.  2-3)  The  Turbulent  Film  Cooling  Boundary  Layer  Region 

The  Investigations  concerning  the  turbulent  boundary  layer  region 
have,  for  the  most  part,  been  experimental  and  have  outnumbered  those 
for  the  potential  core  and  wall  jet  regions.  This  is  primarily  due  to 
the  fact  that  the  turbulent  region  exists  over  the  major  portion  of  the 
area  protected  by  the  film  cooling  blanket.  A  secondary  reason  for  the 
Interest  in  this  region  is  that  here  the  parameters  of  the  flow  will 
most  effect  the  wall  tenperature.  In  the  potential  core  region  sufficient 
coolant  exists  to  prevent  excessive  temperatures  from  being  reached, 
and  If  the  mainstream  velocity  Is  greater  than  the  coolant  velocity, 
the  wall  jet  region  does  not  exist.  Since  this  is  most  often  the  case 
in  combustion  chambers  the  boundary  layer  region  is  of  special  import¬ 
ance. 

The  experimental  investigations  may  be  catagorlzed  into  two  divisions, 
those  with  low  mainstream  temperatures  (<700  °R)  with  film  cooling  or 
heating,  and  those  with  mainstream  temperatures  In  the  range  of  those 
found  in  gas  turbines.  Fortunately,  division  of  the  empirical  results 
into  these  two  groups  eliminates  many  of  the  studies  which  present  cor¬ 
relations  with  little  physical  insight.  However,  one  such  correlation 
of  data  remains  (Papell  and  Trout,  1959)  and  will  be  discussed  first  in 
the  section  below.  The  need  for  a  theoretical  approach  Is  dramatically 
evident  from  such  a  study.  The  work  of  Spalding  (1965)  mentioned  in 
the  previous  section  will  be  briefly  extended  to  the  situation  of  film 
cooling  In  the  turbulent  region.  The  results  of  the  theory  will  be 
compared  to  the  empirical  results  of  Seban  and  Back  (1962).  Finally,  the 
results  of  Hartnett  et  al .  (1961)  will  be  shown  to  be  of  the  same  form 
as  those  predicted  by  Spalding.  The  emphasis  In  this  section  is  again 
placed  on  the  theoretical  approach  since  the  conditions  of  interest 
are  not  adequatly  covered  by  empirical  results. 

Papell  and  Trout  (1959)  were  among  the  first  investigators  to  study 
film  cooling  of  an  adiabatic  flat  plate.  The  range  of  conditions  of 
their  experiments  Is  shown  in  Table  3. 


Table  3. 

Experimental  Conditions 
for  the  Study  of  Papell  and 
Trout  (1959) 


Mainstream  Gas  Temperature 
Mainstream  Mach  Number 
Cooling  Air  Temperature 
Cooling  Air  Mach  Number 
Slot  Height 


520  ®R  to  2000  °R 
0.15  to  0.70 
540  °R  to  870  °R 
0.0  to  1.0 

1/2,  1/4,  1/8,  and  1/16  Inches 
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tflth  very  little  Introduction  as  to  why  tha  efficiency  cf  the  ft  .a 
cooling  slot  should  depend  on  any  specific  quantities,  the  authors  pre¬ 
sented  a  general  correlation  fora: 


where 


Although  the  dependence  of  the  efficiency  on  mass  flow  rates  and  slot 
height  Is  evident  the  authors  provided  no  reasons  for  the  Inclusion  of 
the  temperature  ratio.  Also  absent  from  the  analysis  was  a  discussion 
of  the  physical  reasons  for  organizing  the  components  In  the  manner  shown. 

In  lieu  of  such  empirical  correlations,  a  theoretical  efficiency 
has  been  provided  by  Spalding's  "Unified  Theory".  The  development  of 
a  theoretical  approach  Is  of  Interest  for  two  reasons:  the  results  pro¬ 
vide  Insight  Into  the  phenomem  concerning  the  flow,  and  the  results 
may  more  safely  be  extended  to  conditions  which  have  not  been  experi¬ 
mentally  investigated.  The  relation  between  the  physical  and  mathematical 
findings  have  been  outlined  In  the  previous  section.  The  only  point 
which  should  be  made  Is  that  the  freestream  velocity  Is  no  longer  zero. 
Interaction  of  the  freestream  and  the  boundary  layer  leads  to  the  tur¬ 
bulent  mixing  and  heat  transfer. 

Therefore,  the  results  of  Spalding's  theory  for  film  cooling  effect¬ 
iveness  of  a  turbulent  boundary  layer  may  be  expressed  as: 
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(4-53) 


The  value  of  x  necessary  for  the  Integration  of  Rex  will  not  be  the  dis 


n 


tancfc  from  the  slot  exit,  since  the  boundary  layer  appears  to  start 
somewhere  ahead  of  the  slot  exit. 

The  above  may  be  compared  to  the  experimental  work  of  Seban  and 
Back  (1962),  for  which  Spalding  determined  ■  0.903,  l  •  8.77, 
and  n  *  0.63  for  a  slot  height  of  1/16  Inch.  The  bracketed  term  In 
(4-53)  then  becomes  equal  to  0.210,  and  (4-53)  reduces  to: 


hms  "  hs 
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ps  us  hs 


(4-54) 


The  comparison  between  the  theoretical  and  experimental  efficiencies 
Is  shown  in  Fig.  20  as  a  function  of  (ps  u$/pms  ums)(y|r)S/,52)*  the 
momentum  thl  uness.  The  two  extreme  cases  of  n  9  1 .0  (momentum  transfer 
equals  heat  and  mass  transfer)  and  n  ■  0  (Infinitely  fast  heat  and  mass 
transfer)  are  also  shown  In  Fig.  20. 

A  parallel  analysis  by  Hartnett,  Blrkebak,  and  Eckert  (1961)  pro¬ 
vides  a  similar  relation  to  that  of  Spalding.  These  authors  proposed 
that  the  energy  contained  at  any  boundary  layer  cross  section  of  the 
flow  downstream  of  a  slot  must  be  equal  to  the  energy  added  through  the 
slot, 

^lOT  ■  /oSBL‘)UCPL  <T-T»s>*  (4'55) 

where  L  Is  the  plate  width  and  1SLQT  Is  the  energy  flux.  The  above 
may  be  expanded  to  give: 
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Experimental  results  from  the  work  of  the  authors  showed  that: 
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Similarly 


the  boundary  layer  temperature  profile  was  found  to  be: 


0.2  0.4  0.6  0.8  1.0 
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FIGURE  20  FILM-COOLING  EFFECTIVENESS 
RELATED  TO  MOMENTUM  THICKNESS  (SPALDING,  1968) 
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T  -  T  ^  ,  y  J  3/6  __  _  .  u  2_ 

T - rT —  "  L-C  (-t — )  ]  =  exp  [-0  j 

'aw  'ms  °BL  PBL 

Inserting  (4-57)  and  (4-58)  Into  (4-56)  yields: 


(4-58) 
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Expanding  exp  [-C  (y/<53Lr]  In  a  Taylor  Series  and  retaining  only  the 
first  two  terms: 
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and  note  that. 
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Then 
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However,  the  energy  flux  through  the  slot  may  be  expressed  as: 
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Defining  n  as  before: 
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(4-67) 


The  value  of  K!j  determined  from  experiment  yielded: 

n  *  5.77  (Pr)2/3(Res)-2(^r^-^)“-8  (4-60) 

5  s  ps  us 

which  Is  the  sane  form  as  that  of  Spalding  (Eqn.  (4-54)),  except  for 
the  Inclusion  of  the  Prandtl  number. 

In  conclusion,  the  necessity  for  a  theoretical  method  has  been 
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shown  by  r»v*» winn  the  results  cf  rape'll  and  Trout.  The  cneoretlcal 
analysis  was  the  result  of  Spalding’s  “Unified  Theory"  which  was  sub” 
sisrstlateu  by  the  results  of  Seban  and  Back.  The  Independent  develop¬ 
ment  by  Hartnett,  Blrkebak,  and  Eckert  agrees  with  the  form  of  Spalding's 
result.  Again,  the  emphasis  has  been  on  the  development  of  a  theory 
rather  than  the  use  of  correlations  because  of  the  necessity  of  extending 
the  results  to  higher  temperature  application. 

B.  3)  Penetration  Jets  and  Convective  Cooling 

Although  no  specific  studies  have  been  reviewed  concerning  pene¬ 
trate  t:  a  jets  and  their  effect  on  convective  cooling,  two  authors  (Seban 
et  al.,  1957*  and  Papell,  1960)  have  Investigated  the  efficiency  of 
normal  slot  Injection  and  the  efficiency  of  a  series  of  holes.  Even 
though  the  conditions  covered  In  the  work  of  Seban  et  al.  (1957)  do  not 
correspond  to  those  found  In  typical  gas  turbines,  their  results  Indi¬ 
cate  the  loss  In  cooling  effectiveness  due  to  normal  Injection.  The 
most  ambitious  work  concerning  the  normal  Injection  slot  has  been  done 
by  Papell  (1960),  Included  In  the  study  was  the  effect  of  cooling  with 
a  series  of  small  holes  injecting  normally  Into  the  main  flow.  Due  to 
the  lack  of  Information  on  penetration  jets,  It  Is  believed  that  these 
results  will  indicate  the  trends  of  cooling  efficiency  due  to  normal 
injection  of  air. 

Seban  et  al.  (1957)  studied  the  normal  slot  Injection  into  the 
mainstream  and  Its  effects  on  file  .ooling  efficiency.  The  conditions 
are  similar  to  Wieghardt's  studies  (Stollery  and  El-Ehwany,  1965)  for 
the  heating  of  airplane  wings  to  prevent  Icing.  However,  the  reduction 
of  efficiency  for  csormal  Injection  from  that  of  tangential  Injection 
would  be  expected  to  be  of  the  same  order  as  for  other  mainstream  and 
slot  conditions.  Fig.  21  indicates  a  drop  In  efficiency  of  up  to  59% 
due  to  normal  injection  and  should  Indicate  the  trends  expected  for 
work  at  higher  temperatures. 

Papell  (1960)  investigated  the  effect  of  normal  Injection  In  two 
systems:  the  slot  and  a  series  of  snail  holes  (Fig.  22).  The  conditions 
covered  are  outlined  In  Table  4.  Note  the  fact  that  various  angles  of 
injection  were  used  for  the  slot.  The  angle  v;as  varied  primarily  as  a 
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Table  4. 

Expert  mental  Conditions  for  the 
Study  of  Pape  11  (1960) 


Slot  Configurations 

T-.  *  1500  “R 

01$ 

Ts  -  780  °R 

Hole  Configurations 

Four  rows  Two  rows 
of  1/4  In.  of  1/4  in. 
holes  holes 

Injection 

Angle 

90° 

80° 

45  6 

90° 

90° 

Mai nstream 
Mach  Number 

m 

D 

0.5 

0.3; 

0.6; 

0.7 

0.7 

u 

0.53 

0.75 

0.76 

0.49 

Range 

to 

to 

to 

to 

to 

us 

17.47 

8.38 

9.05 

21.83 

12.80 
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b)  Normal  Hole  Injection  Plates 

FIGURE  22  EXPERIMENTAL  SYSTEM  AND  PLATE  CONFIGURATION  OF  PAPELL  (i960) 
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means  of  confirming  the  physical  model  proposed.  Papell  suggested 
that  the  angle  of  Injection  and  the  final  angle  of  the  flow  or  effective 
Injection  angle  1  could  be  related  utilizing  the  sys  of  vsctsrlal 
velocities  of  the  slot  and  malnstreaai  flow  as  shown  In  Fig.  23.  The 
results  from  such  a  suronatlon  Is  given  by: 
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.  ,  -'us 

c“1+TpTD7 


M-«9) 


In  correlating  the  results,  Papell  found  that  the  cosine  of  the  effective 
angle  gave  the  best  results.  The  effect  was  then  added  to  the  film 
cooling  correlation  used  In  the  previously  cited  work  of  Papell  and  Hatch 
(1961): 

h„,  1  x  K  0.125  IL. 
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(4-70) 
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One  variation  which  should  be  noted  In  Hqn.  (4-70)  Is  the  definition  of 
the  efficiency,  n: 
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FIGURE  23  VECTORIAL  SUMMATION  OF  VELOCITY  COMPONENTS  TO  GIVE 
THE  EFFECTIVE  ANGLE  OF  INJECTION  (PAPELL,  1960) 
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t  =  static  temperature  In  degrees 
T  =  stagnation  temperature  In  degrees 


The  accuracy  of  this  correlation  was  found  to  be  approximate  10*.  Fig. 
24  shows  typical  variations  In  the  efficiency  as  the  slot  angle  was 
varied  with  different  flow  rate  ratios  for  a  1/4  Inch  slot.  These  re¬ 
sults  should  Indicate  the  variation  In  efficiency  due  to  the  normal  in¬ 
jection  rather  than  tangential  injection. 

An  additional  configuration  of  Interest  for  providing  coolant  to 
the  chamber  walls  was  presented  In  the  previously  cited  work  of  Papell 
(1960)  and  will  be  discussed  briefly.  Papell  proposed  the  use  of  a 
series  of  small  holes  as  shown  In  Fig.  22.  As  for  slots,  an  effective 
slot  height  found  empirically  was  used  In  a  correlation  similar  to 
Eqn.  (4-70): 


where 


A$  =  coolant  gas  flow  area 
L|  =  width  of  first  row  of  holes 


For  ums/us  >  1.0  the  correlation  found  was: 


h  LX 

los« n  *  '  cTsrc^7 


0.04]  ( 


^Jk,0-'25  f(k, 

s  s 


Hits 

+  loge  cos  0.8  1ftff  -  0.08  (-~) 


(4-73) 


Fig.  25  Illustrates  typical  results  for  the  above  correlation. 

Typical  of  ths  early  work  done  In  the  field  of  film  cooling,  other 
than  the  model  presented  very  little  physical  Insight  Is  given  as  to 
why  the  correlation  should  be  of  the  form  of  Eqn.  (4-70)  and  (4-73). 
Therefore,  the  above  results  are  to  serve  as  an  example  of  what  the 
trends  are  expected  to  be  for  the  convective  contribution  from  pene- 
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FIGURE  24  EFFECT  OF  VARING  COOLANT  INJECTION  ANGLE  (PAPELL,  1960) 
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FIGURE  25  CORRELATED  RESULTS  USING  EQUATION  (4-73) 
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etratlon  jets. 

5.  4)  Porous  Wails  for  Transpiration  Cooling 

The  concept  of  utilizing  the  effusion  of  cold  annulus  air  through 
a  chamber  wall  made  of  porous  materials  has  been  considered  by  Tlpler 
(1955).  As  with  film  cooling,  the  object  Is  to  provide  a  sufficient 
air  flow  to  form  a  blanket  of  air  to  protect  the  wall  from  convective 
heating.  Although  theoretical  and  empirical  Investigations  have  been 
done  (Eckert  and  Llvlngood,  1953,  Spalding  et  al.,  1964,  Milford  and 
Spiers,  1961,  Spalding,  1965),  they  will  not  be  considered  In  the  light 
of  the  arguments  of  Tlpler  (1955)  and  Milford  and  Spiers  (1961),  which 
are  based  on  four  reasons: 

1)  porous  materials  have  a  higher  exposed  surface  area,  due  to 
the  nature  of  the  material,  which  makes  such  materials  more  prone  to 
oxidation; 

2)  the  porous  material  must  be  cooled  to  a  lower  temperature  than 
solid  walls  for  equal  life,  thus  somewhat  offsetting  the  advantages  of 
having  more  continuous  cooling; 

3)  the  fine  passages  are  prone  to  obstruction  from  impurities  in 
the  Inlet  air  such  as  dust  or  engine  oil; 

4)  the  porosity  of  the  material  may  vary  appreciably  within  a  single 
section  of  wall  leading  to  uneven  cooling  and  excessive  wall  stress. 
Therefore,  the  method  of  cooling  by  effusion  has  been  eliminated  from 

the  present  study. 

C.  External  Convective  Cooling,  C2 

As  for  Internal  flow,  the  external  convective  flow  is  assumed  to 
be  fully  developed  turbulent  flow  and  resemble  flow  through  a  pipe 
(Lefebvre  and  Herbert,  1960,  Winter,  1955,  and  NREC  (Anon.,  1968)). 

Winter  provided  the  following  formulation  of  the  effect  of  external 
flow  with  constant  wall  temperature: 

,  0.214  0.786 

C2  .  0.0255  (^)  (p  Ufl  CP>  <T„  -  TcaI  (4-74) 

where  the  Mussel t  number  Is: 
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and  Is  similar  to  that  found  by  Lyon  and  modified  by  Seban  and  Shlmazakl 
(McAdams,  1954)  for  fully  developed  turbulent  flow  with  constant  wall 
flux: 


hCdA 
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dA  p  u. 

5.0  +  0.025  (~— 
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The  factor  of  5  Is  eliminated  because  it  Is  almost  negligible  compared 
to  the  magnitude  of  the  other  terms.  The  other  two  formulations  use 
studies  done  specifically  In  pipe  flow. 

Lefebvre  and  Herbert  (1960)  and  NREC  (Anon.,  1968)  have  approached 
the  problem  is  the  same  fashion  and  differ  only  slightly  In  the  final 
outcome.  The  studies  of  fully  developed  turbulent  pipe  flow  of  Blalokoz 
and  Saunders  (1956)  were  used  by  Lefebvre  and  Herbert  in  the  following 
form: 
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NREC  (Anon.,  1968)  calculated  C2  with  a  similar  relation  formulated  by 
Humble  et  al .  (1951): 
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Inspection  of  Eqn.  (4-76)  and  (4-77)  reveals  that  they  differ  by  a  factor 
of  (1/Pr)0,4  included  In  Eqn.  (4-77).  For  air  In  the  range  of  tempera¬ 
tures  considered  In  gas  turbines  this  factor  ranges  from  1.10  to  1.15. 

The  relative  ease  of  calculating  either  expression  favors  the  use  of 
the  somewhat  more  precise  value  of  Humble  et  al.  The  assumption  that 
the  flow  Is  fully  developed  and  resembles  the  flow  Inside  of  a  pipe  is 
deemed  a  reasonable  approximation  for  the  calculation  of  Cj. 


CHAPTER  V 

SUMMARY  AND  FUTURE  EFFORTS 

Utilizing  the  superpositioning  theory,  each  of  the  terms  of  Eqn. 
(2-35)  have  been  determined  as  a  function  of  the  wall  temperature. 

R,  +  Cj  +  K  -  R3  +  c2  +  R2  (2-35) 

These  equations  are  summarized  In  Table  5. 

The  solution  of  Eqn.  (2-35)  Is  not  of  supreme  Importance  at  this 
point  in  the  present  study.  Of  more  importance  is  a  summary,  firstly, 
of  what  has  been  accomplished  and  secondly,  what  terms  in  Table  5  need 
to  be  studied  further. 

In  the  first  section  of  the  report,  the  phenomena  present  in  the 
combustion  chamber  were  presented  through  the  use  of  the  conservation 
equations  (Eqn.  (2-1)  -  (2-4)).  A  formal  solution  utilizing  these 
equations  was  shown  to  become  an  Intractable  problem  In  the  presence 
of  the  velocity  and  concentration  profiles  existing  In  the  chamber. 

The  theory  of  radiative  transport  was  then  developed  to  yield  the  ap¬ 
proximation  made  hy  Einstein  (1963).  The  superpositioning  approximation 
was  thus  shown  to  be  the  only  alternative  to  an  exact  solution. 

The  third  chapter  of  the  report  was  concerned  with  the  radiative 
terms  present  In  the  chamber  and  the  analysis  of  each.  The  approximate 
solution  necessary  for  predicting  flame  radiation  on  a  one  and  two 
dimensional  basis  was  reviewed.  The  unfortunate  dependence  of  the  cal¬ 
culation  on  the  empirical  emisslvltles  was  also  revealed.  This  field 
will  most  likely  have  enjoyed  a  great  deal  of  research  since  1965.  One 
objective  proposed  would  be  to  establish  theoretical  reasoning  for 
luminous  emisslvlty  values;  however,  unless  a  scheme  for  predicting  the 
distribution  of  carbon  particles  has  been  provided  In  the  literature 
or  can  be  deduced,  thl s  endeavor  wl  1 1  be  fruitless.  Additional  research 
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Table  5. 

Heat  Transfer  Coop  orients 


R1  "  }  ef  °$B  Tf1,S  (Tf2*5  -  Tw2*5) 
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will  be  conducted  to  determine  tte  advancmaents  and  possibilities  of 
utilizing  tte  theory  of  radiative  transport  In  tte  combustion  chanter. 

In  conclusion,  advancements  In  determining  flans  amlsslvlttes  and  the 
apppllcatlon  of  tte  radiative  transport  theory  appear  to  be  tte  only 
steps  short  of  solving  tte  nore  exact  Integro-differentlal  equation  of 
energy  conservation. 

As  Lefebvre  and  Harbert  (1960)  pointed  out,  tte  convective  contri¬ 
bution  In  the  primary  zone  is  poorly  understood.  Sufficient  knowledge 
of  tte  recirculation  flows  In  tte  prlaary  zone  is  not  available  for  a 
more  realistic  analysis  to  be  considered  (see  Hanond  and  fellor,  1971). 
Hopefully,  tte  lack  of  accuracy  In  detsralnlng  tte  profiles  In  this 
region  will  not  be  teraful  because  tee  annulus  flow  Is  tte  coldest 
around  this  region.  Thus,  the  wall  temperature  rise  should  be  Halted. 
Judging  froa  tte  raaarous  Investigations  of  file  cooling*  this  area 
will  hopefully  provide  useful  advances  especially  concerning  theoretical 
approaches  such  as  Spalding's  (1965).  Finally,  as  aMtlomd  In  tte  sac* 
tlon  on  external  cooling,  the  analogy  between  tte  external  flaw  and  tte 
flow  over  a  flat  plate  Is  believed  to  be  accurate.  Therefore,  If  ad¬ 
vancements  are  to  be  made,  tte  analysis  of  convective  cooling  and  heating 
will  be  tte  most  likely  fields  for  these  advances. 

Although  tte  second  chapter  has  seemingly  ruled  out  the  nore  exact 
approach,  the  approach  of  solving  tte  conservation  equations  will  be 
considered  In  tte  future.  At  this  time,  tte  possible  paths  for  analysis 
appear  to  lead  in  two  directions:  1)  tte  exact  approach  through  solu¬ 
tion  to  the  conservation  equations,  and  2)  approximate  methods  such  as 
the  theory  of  superpositioning.  Tte  literature  review  of  articles 
written  since  1966  should  provide  Intelligent  direction  to  the  most  ac¬ 
curate  approach  possible. 
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